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augment  the  general  description  nf  the  program  provided  in 
the  User's  Manual  and  Mathematical  Models  and  provide  a  pro¬ 
grammer  with  assistance  in  understanding  the  internal  workings 
of  the  SEADYH  programming.  Instructions  for  converting  the 
program  to  various  machines  and  for  modifying  the  program  are 
also  provided. 
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1.0  INTRODUCTION 


The  SEADYN  computer  program  is  a  multifunctional  tool  for  analyzing 
the  structural  response  of  cable,  truss,  and  mooring  systems.  The 
purpose  of  this  manual  is  to  document  some  of  the  major  aspects  of  the 
program  structure.  Program  doctunentation  also  includes  a  theoretical 
toanual  (Ref  1)  and  a  user's  manual  (Kef  2).  The  information  contained 
in  this  manual  is  intended  for  »  computer  programmer  and  the  more 
sophisticated  user  who  may  find  it  necessary  to  get  involved  with  the 
internal  workings  of  the  program. 

The  SEADYN  program  traces  its  pedigree  from  some  early  work  done  by 
Leonard  (Ref  3)  through  a  study  project  carried  out  by  the  Bechtel 
Corporation  for  the  Electronic  Systems  Division  of  the  General  Electric 
Company  (GE).  The  end  product  of  that  work  was  a  program  called  NLIN 
(Ref  4).  This  author  was  technical  advisor  on  that  project.  The  SEADYN 
program  began  where  NLIN  left  off  as  part  of  the  author's  doctoral 
project  in  1974.  It  has  been  under  various  levels  of  development  since 
then  with  sponsorship  coming  from  GE,  the  Chesapeake  Division  of  the 
Naval  Facilitie.s  Engineering  Command  (CHESDIV),  and  finally  the  Naval 
Civil  Engineering  Laboratory  (NCEL) .  The  support  of  CHESDIV  led  to  the 
implementation  of  the  ship  mooring  capabilities  for  static  and  frequency 
domain  solutions  in  1976.  NCEL  support  has  provided  for  the  pay-out/reel-in 
capability,  the  drag  amplification  calculations  for  strumming,  the 
viscous  relaxation  sialic  solution,  the  bottom  limited  catenary  element, 
variou.^  revisions  and  improvements  in  the  solution  algoritlims,  and  a 
complete  restructuring  of  the  program  including  the  new  free-field  input 
format . 

A  plotting  po.st-processor  program  which  is  based  on  the  restart 
file  structure  described  in  this  document  is  also  available  (!>]. 
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2.0  PROGRAM  STRUCTURE 


T){£  SEADYM  program  is  composed  of  a  simple  main  program,  which 
establishes  the  size  of  the  working  storage  in  COMMON/ACOM/  and  numerous 
subroutines.  The  heart  of  the  program  is  found  in  the  SEADYN  and  MANIPR 
subroutines.  SEADYN  is  called  by  the  main  program  and  manages  the  setup 
of  storage  partitioning,  reading,  and  deciphering  of  the  free-*field 
input  and  the  initial  access  to  the  restart  files.  Control  is  then 
passed  to  the  MANIPR  subroutine  to  manage  the  actual  problem  solution 
through  calls  to  the  subanalysis  routines. 

The  programming  effort  has  attempted  to  maintain  modularity  in  the 
structure.  Major  subanalyses  are  handled  in  separate  groups  of  subroutines 
that  are  called  in  response  to  input  requests.  The  overall  analysis  and 
the  individual  solution  options  have  been  written  to  be  nearly  independent 
of  the  type  of  element  used  to  generate  the  equations.  For  example,  the 
clement  stiffnesses  and  internal  force  components  are  generated  from 
calls  to  a  single  controlling  routine  regardless  of  the  type  of  analysis 
being  performed.  This  modularity  should  prove  to  be  useful  if  the 
addition  of  new  subanalyses,  solution  options,  or  element  types  are 
contemplated. 

The  program  is  written  entirely  in  FORTRAN  IV  (CDC-Extended  FORTRAN). 
Initial  developments  were  on  a  GE-635  computer  with  compatibility  with  a 
CDC-6600  computer  maintained.  Since  1977,  development  has  been  on  a 
CDC-7600.  Some  effort  has  been  directed  at  restructuring  to  facilitate 
conversion  to  ANSI  standard  FORTRAN-77  and  to  other  machines.  Conversions 
of  interim  versions  have  been  accomplished  on  the  VAXll/780,  PRIME,  and 
DATA  GENERAL  machines.  The  major  items  encountered  in  conversion  include: 

1.  The  free-field  reader  —  uses  ENCODE/DECODE 

uses  REFORMAT  character  manipulations 

2.  The  CDC-prograra  statement  in  the  main  program 

3.  End  of  file  checks  using  IF(EOF)  a,b 

A.  The  use  of  character  strings  in  titles,  date  and  time  (including 
Hollerith  word  lengths) 

‘S.  Asterisk  (*)  comment  cards 

ti.  .Asterisk  (*)  and  quote  (”)  format  delimitors 

7.  The  need  for  double  precision  on  32-bit  word  ’machines 

8.  .Seven  character  subroutine  name:  RESTART 

The  global  logic  structure  is  represented  in  Figure  2.1.  The 
relationships  of  the  major  subroutines  are  outlined  in  Figure  2.2. 
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Figure  2.1.  Macro  flow  chart  for  SEADYN. 
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Figure  2.2  Outline  of  major  subroutine  relationships. 


3.0  DESCRIPTIONS  OF  COMMON  STORAGE 


Common  storage  can  be  segregated  into  four  categories;  working 
area  for  variable  dimension  data,  fixed  length  storage  for  restart  data, 
scratch  storage  for  communication  between  subroutines,  and  special  data 
tables.  Each  of  these  are  described  below. 


3.1  Working  Storage  for  Variable  Dimension  Data 

COMMON/ ACOM/  A(XXm) 

The  length  of  the  A  array  is  specified  in  the  MAIN  program.  Partition¬ 
ing  of  A  into  the  data  arrays  needed  In  the  solution  is  done  by  the 
SEADYN  and  h'ANIPR  subroutines.  SEADYN  reads  the  problem  description 
data  or  restart  file  to  determine  the  size  of  the  data  arrays.  SEADYN 
then  passes  the  first  position  of  each  array  to  MANIPR  through  its 
calling  sequence.  With  only  two  exceptions,  /ACOM/  is  referred  to  only 
indirectly  through  the  subroutine  calling  sequences.  The  two  exceptions 
are  the  MODES  and  FREQ  subanalyses  which  place  the  contents  of  /ACOM/  on 
a  scratch  file  using  RESTART  and  then  use  the  storage  for  other  purposes. 
/ACOM/  is  restored  through  RESTART  upon  exit  from  these  subanalyses. 

The  partitioning  of  /ACOM/  is  described  in  Table  3.1. 

3.2  Fixed  Length  Storage  for  Restart  Data 

In  addition  to  the  element  and  node  data  contained  in  /ACOM/,  the 
essential  problem  data  are  contained  in  a  set  of  common  blocks  with 
fixed  dimensions.  These  are  arranged  in  pairs  which  contain  floating 
point  and  fixed  point  data  in  segregated  forms.  A  separate  common  block 
contains  all  of  the  logical  variables.  These  blocks  are  written  and 
read  for  restart  purposes.  The  segregation  is  made  to  facilitate  the 
conversion  to  machines  which  use  byte-oriented  storage  with  type-dependent 
word  structures. 

The  use  of  the  common  blocks  is  summarized  below.  The  actual 
structure  of  the  common  blocks  is  given  in  Table  3.2. 


/BUOYS/,  /IBUOYS/ 

/CABLE/,  /ICABLE/ 

/CONTRL/,  /ICNTRL/ 
/DSPCON/,  /IDSPCN/ 

/PAYOUT/,  /IPAYOT/ 
/SHIPS/,  /ISHIPS/ 
/STRUM/,  /ISTRUM/ 
/TIMED/,  /ITIMED/ 
/LOGIC/ 


Body  tables,  limit  set  data,  body  and  limit 
location  data,  surface  buoy  data  and  related 
control  data 

Line  element  data  and  material  tables,  fluid 
data 

Basic  control  information  for  the  solution  options 

Static  imposed  displacement  data  and  norm  data 
for  static  iterative  solutions 

Pay-out/ reel -in  data 

Data  for  ships  and  static  ship  loads 

Strum  string  data 

Transient  response  data 

Logical  variables  used  for  solution  control 
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Table  3.2.  Fixed  Length  Restart  Coaaon  Blocks 


COMMON/BUOYS/ 

1  AABUOY(IOO) ,ADM(50) ,BAMC(50) ,BLEN(50) ,BMOM(50) ,B0MAS(50) . 

2  BOV0L(50),BSCD(50),BU0KP(3,50),BWND(50),CBOO,CORDLM(50), 

3  CYLRMS(50) ,DBU(50) .FBUSR(50) ,RELFAC(50) ,SBAMP(50) »SGHDLD, 

4  TOLIM(50),UBS(50),VIB(3) 

COMMON/ I BUOYS/  lAABU(lO). 

1  IBS(50) ,IBU(50) ,IDRB(50) ,IMPB0D,IMPN0D,I0PT,JANCR(50) , 

2  JSLP (2 , 50) .KNSTRN ,KONECT( 10 ,50) , LB0DN(50) ,LIMN0D(50) , 

3  LIMSET(50) ,MBLN(50) ,MEDMB(50) ,HORBUO,MXBLOC,HXBODY, 

4  MXLIMS,MXLLOC,NBLOC,NBUTAB,NCYLB,NLIMS,NIXOC,NMOTN(50) 

COMMON/CABLE/ 

1  AACAB(100),CABMAS(10),CAMC(10),CURMUL,DIAM(10),E(20,10), 

2  FDEPTH(5),FGAM(5),FLPAR(10,10),FVISC(5),G3(10),STR(20,10), 

3  TENULT(10),TT(20,10),TTD(10),TTK(10) 

COMMON/ICABLE/  lAACA(lO). 

1  IDAMPR,IDRG(10),IFLCOD(10),MATDMP(10),MATT,ME(10),MED(10), 

2  MXFLOW,NFLUr,NFLVRY 

COMMON/ CONTRL/ 

.  AACON(IOO) ,ACCFAC,DELFAC, DERAD, DELTMP, DERR, DMU, 

1  EXTRAP , FDI VY , FINC , FLNVY , FRCVY , GRAY , G1 , G2 , OVRSHT , 

2  PARMT , PI , PINC , PSTEP ,RATD ,RERR , SRCHFC , SSTART 

COMMON/ I CNTRL/  lAACL(lO), 

3  I BEND , I BEND 1 , I BEND2 , IBG , IDIR , IFCNT , IFXFL , IGK , IKNSTN , INC , 

4  INDRAG, INVY,IPR,ISTART,IUPDT,JDLD,JDYN,JMPDT,JOVR, 

5  KONVRT , KOUNT ,KUP , LMITER ,MOSTAT ,MVB ,MVBINC ,NBASE ,NCONC , 

6  NE,NFN,NF3,NFLUID,NFIX,NN,NN3,NPRST,NRUP,NSLAVE, 

7  NSTUP , NTYPE , NUMSET , NUP , NXTYPE , MODEI 1 , M0DEI2 

COMMON/DSPCON/ 

1  AADSP(IOO) ,DISPC(30) ,DISPP(30) ,HIRSDL,RNORM,RNORMP,RNRMPP, 

2  VNORMP,\’NPP,VNPPP 

COMMON/ IDSPCN/  lAADS(lO), 

1  IDIS(30,3),IDWN,TNUP,ITUP,JSTEPR,LMKEEP,NDISP,NDR0P, 

2  NFKEEP(30),N0DWN,NSL0P 

COMMON/PAYOUT/ 

1  AAPAY(IOO) , AMAXL(5) ,CURLEN(5) ,DELT,DLREr(5) ,DPTMF(5) , 

2  ESP(2,5),PAYV(5),PTMF(5),TFSAV(5),UMVB(15,3),UMVBP(15,3) 

COMMON/ IPAYOT/  lAAPO(lO), 

1  J0P(5) , JPELT(5) ,MITN0T,MULMAT(5) ,MULTIM,NELP0I(5) ,NGR0W(5) , 

2  NNPOI(5),NOP,NPOVRY(5),NSHRNK(5) 
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Table  3.2.  (Continued) 


COMMON/SHIPS/ 

1  AASHP (1 00 ) , ACCCUR , ACCWND , APROP ( 3 ) , BLOCK ( 3 ) , CAD , CPROP ( 3 ) , 

2  CR(3),CS(3),CURCOE(20,3,5),CURHF.D(20),CURin',CURVF,L(5),DRAD, 

3  FACH(3) ,FSFRC(3) ,FSFRW(3) ,FSLEN(3) ,FSVEL(3) ,GAIR,HEAD(3) , 
HEDEND,HEDIMX,HEDINC,HEDNOW,PROPF(3) ,RATL(3) ,SAE(3> ,SAS(3) , 

5  SBEAM(3) ,SDRFT(3) ,SDSPV(3) ,SFACW(3,3) ,SHIPK(4,3) ,SHPKP(3,3) , 

6  SHTRN(3,3,3),SLT(3),SLWL(3),SURFCE,TSAPL(3),VAIR,WAD,WDEPTH(3), 

7  WI^«),WNDCOE(20,3,b),WNDHED(20),Va«)VEL(5) 

COMMON/ I SHI PS/  lAASHP(lO), 

1  ICR,IFREQ,ISHIP(3),ISHTAP,ISURLD,IUP,KODEC(2),KODEW(2), 

2  LSHP ( 3 ) .MXSHIP ,  .CRNT , NSFILE ,NSHIPS , NTHETC ,NTHETW , NUMHED  ^NWIND 

COMMON/ STP.UM/ 

1  AASTM(100,CEPS.CRST(30),GKS(3,3),QST(20),RVELN,RVELNP, 

2  STLEN(20) 

COMMON/ ISTRUM/  lAASTM(lO), 

1  1 STRNG ( 30 ) , I STROP , KSTRNG (20,30), MXSELT , MXSTRG , NSTRNG , NSTRUP 
COMMON/TIMED/ 

1  AA  T I M ( 1 00 ) , ALPNEW , BETNEW , DALPHA , DBETA , DMAXAB , DMAXP , 

2  DT , UTH , DTL , DTLL , DTMAX , DTRSRT , DTU , ERR , FTF ( 3 ) , FTI ( 1 5 ) , GAMNEW , 

3  T ,  TM.AX ,  TMFFF (  3 ) ,  TMFRM (15), TPAR1>  ( 20 , 20 )  ,  TRSRT ,  TZ ,  iJB  ( 1 5 ) 

COMMON'/ ITIMLO/  lAATIM(lO), 

1  [ALTR,JB(5},TLF(3),IMTMF(15),1MX,ISIGN1,ISIGN2,ITFCOD(20), 

2  I TOP , JB ( 1 5 ) , KNTRST , MXTFUN , N I XPRN 


3.3  ScratcL  Storage  Common  Blocks 

Temporary  data  or  data  not  required  by  restart  are  handled  in  the 
following  common  blocks.  Table  3.3  lists  the  contents. 

/CDCAL/  Temporary  data  for  strum  string  calculations 

/CHKDAT/  Component  check  data  for  element  fluid  loads 

/FRQDAT/  Frequency  domain  interface  data 

/HEDCHR/  Titles  and  labels  for  page  headings 

/HEDDAT/  Control  data  for  page  headings  and  working  storage  for 
free-field  input  processing 

/RETAPE/  Rest.'trt  file  data 

/SHPLBL/  Contains  the  ship  load  title  block  and  labels  for  the 
ship  load  data  set  in  current  use 

/SHP'fAP/  Contains  the  ship  dimension  data  for  the  ship  load  data 
set  in  current  use 

3 . 4  Conuiion  Blocks  for  Data  Tables 

Special  purpose  data  tables  are  contained  in  the  following  common 
blocks.  Table  3.4  lists  the  contents. 

/COMPNT/  Component  inventory  tables 

/SIZE/  List  of  common  sizes 

/T.APF.S/  List  of  file  names 

/TYPES/  List  of  numeric  codes  for  input  keywords 
/IfLAC/  Array  of  input  keywords 

The  mooring  component  inventories  are  contained  in  arrays  which  are 
gene) aiecl  by  data  statements  in  the  COMPDT  Subroutine.  An  explanation 
of  the  entries  in  the  component  inventory  common  block,  /COMPNT/,  is 
given  below. 

DKSCRIPTION 

.ATYPE  Labe]  for  each  of  6  anchor  types.  Each  label  is 

three  6-character  Hollerith  words. 

.ANCTAB  Storage  for  anchor  data  for  up  to  six  anchor  types 

(third  dimension).  The  array  allows  up  to  16  anchors 
of  each  type,  listed  in  order  of  increasing  weight. 
Five  items  are  given  for  each  anchor. 
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NANCK 

BTYPE 

BUOTAB 


NBIIOY 

HTYfE 

BAWTAB 

HAWSIZ 


NllAWS 

CHAIN 

NChN 

HLDFAC 


Word  Contents 

1  Anchor  weight 

2-4  Hollerith  labels  for  federal  stock 

number 

5  Holding  power 

The  holding  power  is  calculated  from  the  weight  and 
the  factors  given  in  HLDFAC. 

An  array  listing  the  number  of  anchors  of  each  type 
in  the  inventory. 

Label  for  each  of  two  buoy  types.  Each  label  is 
three  6-character  Hollerith  words. 

Storage  for  buoy  data  for  up  to  two  buoy  typer,.  The 
array  allows  up  to  six  buoys  of  each  type,  listed  in 
order  of  increasing  buoyancy.  Seven  items  are  given 
for  each  buoy: 

Word  Contents 

Outside  diameter 
Height 
Weight 

Nominal  buoyancy 
Maximum  buoyancy 

Ho’ ;  rith  label  for  federal  stock 
number 

An  array  list  .  he  number  of  buoys  of  each  type 
in  the  inven' 

Label  fo  -c ;  •  four  hawser  types.  Each  label  is 

three  6->.  •  •  >llerith  words. 

Stora*'  sr  data  for  up  to  four  liaswer  types. 

T''''  .  i  up  to  23  sizes  for  each  type.  The 

tWk  ..  ,  -r.  for  each  size  are  the  tensile  strength 

and  the  weight  per  100  units  of  length. 

A  list  of  the  haswer  sizes  given  in  increasing  order. 
It  is  assiuned  that  all  four  hawser  types  have  the 
same  list  of  sizes.  If  no  entry  is  available  for 
a  given  size,  then  set  the  strength  to  zero  and  the 
weight  to  some  small,  nonzero  nujr<ber. 

The  number  of  howser  sizes  for  each  howser  type.  (Set 
to  23  for  each  hawser  type.) 

Storage  array  for  chain  data.  Presumes  oni^  one  chain 
type  (stud-link  chain).  The  array  contains  31  entries 
for  size,  strength,  and  ’’eight  per  unit  length. 

The  number  of  chain  •  les  (31). 

A  list  of  the  holdi.;  power  factors  for  the  six 
anchor  types.  This  is  the  number  which  multiplies 
the  anchor  weight  to  get  the  hoi.  ng  power  in  firm 
sand . 
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The  component  data  presently  listed  in  the  inventories  assume 
weights,  buoyancies,  and  strengths  are  in  units  of  pounds.  Lengths  and 
buoy  dimensions  are  in  feet.  Hawser  and  chain  sizes  are  in  inches. 

It  would  be  possible  to  alter  the  contents  of  the  inventories  by 
rewriting  the  COMPDT  subroutine.  The  subroutine  is  called  once  each 
time  a  component  check  or  design  selection  is  made.  It  is  presently  set 
up  to  avoid  re-calculating  entries  after  the  first  call.  The  program 
uses  the  items  in  the  arrays  to  search  for  the  appropriate  entries.  The 
appropriate  values  are  then  selected,  scaled  without  altering  the  arrays, 
and  used.  Minor  changes  could  be  made  without  changing  the  calling 
program.  More  extensive  changes  or  generalizations  may  require  some 
modifications  (not  major)  in  the  calling  program. 

A  listing  of  the  inventories  is  presented  in  Appendix  D  of  the 
User's  Manual. 
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Table  3.3.  Scratch  Storage  Comoon  Blocks 


COMMON/CDCAL/  AMAX ,OMS ,OMG ,RTIN,YSAVE, 1ST 
COMMON/ CHKDAT/  NELCK,FELT(3 ,2)  .INVFLG 

COMMON/FRQDAT/  SPECA,SPECB,DOMG,OMGMN,OMGMX,AMPMN,FRCFAC,FACLF,N , 

1  TIMFAC ,  IBFG ,  lUNRES ,  IDRITR ,  ICONMS ,  ICMCHF ,  IFRQUP ,  IROLIT 

2  , KF I LEF , NOB ,  NOH , NOK , NRV , MNOB ,MNOH , MNOK , MNRV , ISHPFL , ITCONF , NSOLN , 

3  ELL,TMAS,BAM(30),BA4S(8),DA(6,6),DB(6,6),DC(6,6),GMU(6,6), 

A  HDGl (30)  ,RAMG(8) ,WAVEL(30) ,OMG,NFR£QS,NEWCON,N£WR£D,VDUM(3) , 

5  DRFTFR(3) ,TEMP(6,6) ,GHED,NSRDC,OM2,WVLN,WVAMP»WVSLP,NS3 

COMMON/HEDCHR/  IDAY,ITIME,HED(20) 

COMMON/HEDDAT/  IPAGE,NLINES,MXLINE,NW0R])S,W0RI>S(100)  »NHED, NOLINE 
COMMON/RETAPE/ 1  RST,NTAPE,NFLILE 

COMMON/SHPTAP/  SCALE , TDEPTH .TBLOCK ,TSLT , TSAE , TSAS , TSWL , TSB , TSK , 

1  TSDSP.TSAP 

COMMON/  SHPLBL/  SHPCAP (12), WLBL , CLBL , LLBL , VLBL 
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Table  3. A.  Common  Blocks  for  Data  Tables 


COMMON / COMPNT/  ATYPE ( 3 , 6 ) , ANCTAB (5,l6,6),NANCR(6),DTYPE(3,2), 

1  BUOTAB(7 ,6,2) ,NBU0Y(2) ,HTYPE(3,4) ,HAWTAB(2,23 ,A) ,HAWSIZ(23) , 

2  NHAWS(4),CHAIN(31,3),NCHM,KLDFAC(6) 

COMMON/TAPES/  NIN.NOUT.NTAPEl ,NTAPE2 
.  ,NTAPE9,NTPCHK 

COMMON/SIZE/  NINA, NINE, NINRL,N1NDSP,NINP0,NINSIIP,NINSTM,NINT, 

1  NINC ,NCOM, IFILE(4) ,NPRECZ,N1NBI ,NINCI .NINCLI ,NINDSI , 

2  NINPOI ,NINSHI .NINSTI .NINTMI 

COMMON/TYPES/ 

1  DE AD , AL I VE , DYNAM , TSSS , AMODE , FREQ , CHEK , PLOT , ANEW , F IN I S , 

2  PROB , REST .BLOC , BODY ,ELEM , FLOW, FlUI ,AINVE ,ALIMI , ALINE , 

3  ALLOC , AMATE , ANODE , ASHI P , ASTRUM , TENS , TFUN , CURR , FIX , FREE , 

4  AIMPA,AINIT, IKEEP.ALOAD, ALVAR, AMOV'E.AMSOL.OUTP.PAYO, SAVE, 

5  SBUO , SOLD , ASTEP ,SURF , ATIME , AWIND ,FSOL , SPEC ,EXTE ,RESU , 

6  AHEAD , RAND , REGU , DONF. , ANCH , ABUOY , ACONF 

COMMON/TFLAG/  IFLAGS(60) 
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4.0  STORAGE  REQUIREMENTS  AND  VARIABLE  DIMENSION 

The  amount  of  computer  storage  required  to  successfully  execute  an 
analysis  case  can  vary  with  the  number  of  nodes  and  elestents  as  well  as 
with  the  subanalysis  options  selected.  Direct  restrictions  on  the 
number  of  nodes  and  elements  are  avoided  by  the  use  of  the  implied 
dimensioning  (partitioning)  of  COMMON/ACOM/ .  Since  the  present  version 
of  SEADYN  uses  an  in-core  solver  for  both  the  stiffness  equations  and 
the  eigenvalue  problem,  the  maximum  problem  size  is  restricted  by  the 
amount  of  storage  available.  The  working  space  for  these  solutions  is 
the  /ACOM/  region.  The  stiffness  equation  solution  uses  the  portion  of 
/ACOM/  left  after  the  node  and  element  data.  The  eigenvalue  solution 
copies  /ACOM/  to  a  scratch  file  and  uses  the  full  length  of  /ACOM/  (see 
Section  3.1  and  Table  3.1). 

The  size  of  /ACOM/  is  specified  in  the  main  program,  which  is 
simply  a  starter  deck  to  assign  a  length  to  /ACOM/  and  call  SEADYN. 

Besides  the  CDC  PROGRAM  statement,  the  main  program  contains: 

COMMON/ACOM/  A(xxxxx) 

NC0M=XXm 

CALL  SEADYN  (NCOM) 

STOP 

END 

The  minimum  required  value  for  NCOM  depends  on  the  problem  size, 
and  the  subanalysis  options  selected.  The  amount  of  data  needed  to 
store  the  element  and  node  data  is  computed  by: 

NBASEs:31*NE  +  66*NN 

The  global  stiffness  matrix  (GK)  begins  its  storage  at  NBASE+l.  The 
stiffness  matrix  storage  format  is  by  rows,  starting  at  the  diagonal 
elem.^nt  and  going  to  the  ha  If -bandwidth'-^  (IBEND).  Thus,  any  solution 
which  uces  the  global  stiffness  matrix  will  require  at  least  NBASE  ♦ 

NL'3''IREND  storage  locations  in  /ACOM/  (NF3  was  defined  in  Table  3.1). 

This  amount  is  required  in  the  MNR,  RFB,  VRS/VRR,  and  SLl  solution 
methods.  The  DYN  SAO  with  the  DIM  option  does  not  use  the  global 
stiffness  matrix.  (DYN-SLI  with  moved  nodes  imposes  some  special 
restrictions  since  more  storage  is  needed  for  GK,  and  the  moved  nodes 
must  be  numbered  just  before  the  slave  nodes.  Use  of  this  option  is 
highly  discouraged.  It  is  complicated  to  use  correctly  and  much  less 
efficient  than  the  DIM  solution.) 

The  eigenvalue  solution  (MODES)  uses  the  Jacobi  method  and  starts 
its  working  storage  at  the  beginning  of  /ACOM/.  It  begins  with  a  diagonal 
mass  matrix  (GMM)  and  a  compacted  version  of  the  stiffness  matrix  (GKK) 
which  has  the  full  upper  triangle  of  the  stiffness  matrix.  GKK  requires 
H(NF3)  (NF3-*-!)  words.  Next  is  a  working  area  which  is  NF3*NF3  words.  The 
eigenvectors  (NF3  of  them)  end  up  here  as  the  TGK  array.  The  eigenvalues 
follow  in  the  next  NFS  locations.  The  required  space  is  then  . 5*NF3*(3*NF3+5) 


’^he  half-bandwidth  is  the  largest  number  computed  for  any  row  by  subtracting 
the  row  niimbcr  from  the  number  of  the  rightmost  column  which  contains  a 
nonzero  entry  plus  one. 
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The  frequency  domain  solution  allows  for  the  formation  of  a  consistent 
mass  matrix  (CGM)  which  is  created  immediately  following  GK  in  /ACOM/ . 

Since  it  is  the  same  size  as  GK,  the  storage  needed  is  NBASE+2*NF3*IBEND. 
The  FREQ  solution  uses  complex  arithmetic.  The  system  coefficient 
matrix  is  generated  by  copying  GK  and  CGM  into  a  prior  portion  of  /ACOM/ 
in  complex  format  (see  Table  3.1). 

The  variable  storage  space  needed  for  NCOM  is  summarized  below  for 
each  SAO  type; 

DEAD,  LIVE  TSSS,  DYN 

NCOM  =  NBASE  +  NF3*1B 

where  IB  =  IBEND  for  all  but  DYN-DIM 
IB  =  0  for  DYN-DIM 


MODE 

NCOM  -  .5*NF3*(3*NF3+5) 


FREQ 

NCOM  =  NBASE  +  2*NF3*IBEND 


NCOM  =  NBASE 

Storage  size  checks  are  made  at  the  beginning  of  each  SAO  to  determine 
It  enough  space  is  available.  If  not,  a  message  is  printed  to  indicate 
the  sp.icc  needed,  and  the  run  is  aborted.  Note  that  for  problems  using 
more  than  one  SAO,  the  largest  NCOM  must  be  used. 

The  CDC  version  of  SEADYN  uses  single  precision  logic  except  for 
.iccumulators  in  the  simultaneous  equation  solvers  in  the  SLVBAN  and 
rOMBAN  subroutines  (see  Section  7.0).  On  machines  with  word  lengths 
less  than  60  bits  the  program  should  u.se  double  precision  on  most  of  the 
i loatiug  point  computations.  The  critical  computations  are  the  length 
ihauges,  stiffness  matrix  coefficients,  and  components  of  the  force 
residual.  The  most  direct  .and  reliable  approach  to  treating  double 
nrec  isioii  is  to  type  all  floating  point  variables  double  precision. 

This  is  most  readily  done  with  FORTRAN-77  or  IBM-FORTRAN.  If  this 
approach  is  taken,  then  the  variable  NPRECZ  in  the  SEADYN  subroutine 
should  be  set  to  2.  This  establishes  the  appropriate  niuober  of  words 
to  be  written  and  read  for  the  restart  files.  When  NPRECZ=2  the  vari¬ 
ables  in  the  COMMON  statements  of  the  RESTART  subroutine  should  not  be 
typed  double  precision.  > 

If  it  is  decided  to  selectively  type  the  variables  double  precision 
some  customizing  of  the  variable  dimension  process  will  be  required. 

Tlie  critical  variables  (see  Section  8.0)  are  listed  below: 

ACOM-Cross  Reference  (Table  3.1) 

DS,  DSO,  DSR,  XC,  XO,  XS,  U,  UP,  US,  UD,  UDP,  UDD,  UDDP,  UDDS, 

F2,  GK,  DU 


15 


Spe c lal  Temporary  Variables/Subroutine ; 

CL,  CDX/RESIDL  CL,  CDX/UPDATF. 

Since  the  ACOM  cross  reference  requires  the  A  array  partitioning  to 
be  expended,  the  variables  NBASE  and  NINA  (SEADYN  and  MANIPR  subroutines) 
must  be  increased.  In  addition  it  will  be  necessary  to  modify  the  NC/ . . 
variables  in  the  SEADYN  Subroutine.  The  least  troublesome  approach  (bul 
storage  extravagant)  is  to  make  A  double  precision  and  to  double  NINA 
and  NBASE.  Then  the  NC..  variables  do  not  need  to  be  adjusted.  If  A 
is  typed  double  precisioA  then  NCOM  should  be  doubled  since  it  (like 
NINA  and  NBASE)  count  single  precision  floating  point  words. 

The  arrays  which  do  not  depend  on  the  number  of  nodes  or  elements 
have  been  given  fixed  dimensions.  Should  these  need  to  be  changed,  it 
will  be  necessary  to  make  three  types  of  changes.  They  are: 

1.  Change  the  dimension  of  the  specific  arrays  involved.  A  cross 
reference  between  the  COMMON  involved  and  the  descriptions  of 
the  variables  in  Section  8.0  will  clarify  which  ones  need  to 
bo  changed  (see  Table  4.1). 

2.  Change  the  parameter  which  specifies  the  maximum  number  allowed 
(e.g.,  MXBODY). 

3.  Change  the  parameter  which  specifies  the  total  number  of  words 
in  the  COMMON  block  (e.g.,  NINB) .  This  is  used  by  RESTART  to 
determine  how  many  words  to  read  or  write. 

Care  will  need  to  be  taken  to  assure  the  COMMON  statements  in  all  of  the 
affected  subroutines  are  changed.  It  should  be  noted  that  making  such 
change.*:  causes  the  restart  files  to  be  unique  and  not  accessible  to 
other  version.s  of  the  program  or  to  the  compatible  SEAPLT  graphics 
progiam. 

Table  4.1  is  a  list  of  the  array  dimension  limits  with  the  arrays 
and  paiametors  involved.  Table  4.2  lists  othi;r  program  size  restrictions. 
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Tabic  4.1.  Fixed  Dimension  Array  Data 


Item 

Limit 

Variable 

Name 

Common 

Size 

Arrays 

Involved 

Bodies  in  BODY  table 

50 

MXBODY 

NINB 

NINBI 

BAMC(50) 

BLEN(50) 

BM0M(50) 

BOMAS(50) 

BOV0L(50) 

BSCD(50) 

BWND(50) 

DBU(50) 

CYLRMS(50) 

IDRB(50) 

Body  locations 

50 

MXBLOC 

NINB 

NINBI 

FBUSR(50) 

SBA11P(50) 

UBS (50) 

IBS(50) 

IBU(50) 

JSLP(2,50) 

LB0DN(50) 

MBLN(50) 

NMOTN(50) 

BU0KP(3,50) 

Limit  Conditions  table 

50 

MXLIMS 

NINB 

NINBI 

CORDLM(50) 

RELFAC(50) 

TOLIM(50) 

JANCR(50) 

Limit  Loc.^tlons 

1 

i 

50 

MXLLOC 

NINBI 

K0NECT(10,50) 

LIMNOD(50) 

LIMSET(50) 

1  Li  tun;  ccnnected  to  a  limited 
node 

10 

NINBI 

KONF,CT(10,50) 

Cable  mal.t‘rial  tables 

i 

1 

i 

1 

1 

1 

10 

NINC 

NINCl 

CABMAS(IO) 

CAMC(IO) 

DIAM(IO) 

E(20,10) 

STR(20,10) 

TT(20,10) 

G3(10) 

TENULT(IO) 

TTD(IO) 

TTK(IO) 

IDRG(IO) 

JFLCOD(IO) 

MATDMP(IO) 

ME(IO) 

MED(IO) 
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Table  4.1.  (Continued) 


Limit 


Variable  Common 
Name  Sixe 


Entries  in  any  material  table 


Entries  in  FLUID  table 
(program  logic  also  involved 
here) 

Ship/Platfonn  rigid  bodies 
(program  logic  Units  FREQ 
solution  to  only  one  ship) 


Arrays 

Involved 


E(20,10) 

STR(20»10) 

TT(20,10) 

FDEPTH(5) 

FGAM(5) 

FVISC(5) 


NINSHP  APROPO) 

NINSHI  BL0CK(3) 
CPR0P(3) 

CR(3) 

CS(3) 

CURCOE(20,3,5) 
WNDCOE(20,3,5) 
FACH(3) 
FSFRCO) 
FSFRW(3) 
FSLEN(3) 
FSVEL(3) 
HEAD(3) 
lSHir(3) 
LSHP(3) 
PR0PF(3) 
RATL(3) 

SAE(3) 

SAS(3) 

SBEAM(3) 
SDRFT(3) 
S0SPV(3) 
SFACW(3) 
SHPKP(3,3) 
SI.T(3) 

SI.WL(3) 
TSAP(3) 
WDEPTH(3) 


reel-in  ends  (must  be 
Toss  than  or  equal  to  the 
iiunber  of  moved  nodes  allowed) 


NINPO 

NINPOI 


AMAXL(5) 

CURIJEN(5) 

DLR£F(5) 

DPTrtF(5) 

ESP(2,5) 

PAYV(5) 

PTMF(5) 

TFSAV(5) 

J0P(5) 

JPELT(5; 

MULMAT(5) 

NELPOKS) 


i  >*-w. 


ttrjTTWtyion.-n'.i*  ^a’B  vi 
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Variable  Connoii 
Name  Size 


Arrays 

Involved 


Moved  nodes  in  DYN  or  TSSS 
(Note:  the  row  dimension  is 
3  times  the  number) 


Nodal  components  with  imposed 
displacements  in  DEAD,  LIVE 
or  MODE 


St  nun  strings 

Elements  in  a  strum  string 

Number  of  flow  fields  in  FLOW 

I’.ii  auieV  04  s  p«*r  flow  field 
i  ir,i.-  functions  defined  by  TFUN 

Parameters  per  time  function 

Load  variation  sets  defined 
by  I.OAD/LVAR 


NGK0W(5) 

NNP0I(5) 

NP0VRy(5) 

MSHRNKCS) 

f) 

NINPO 

NINT 

NINTMI 

UMVB(15,3) 

1)MVBP(15,3) 

FTI(15) 

TMFRM(15) 

UB(15) 

IB(5) 

1MTMK15) 

JB(15) 

30 

NINDSP 

NINDSl 

DSPC(30) 

DISPP(30) 

IDIS(30,3) 

NFKEEPOO) 

30 

NINSTl 

1STRN0(30) 

KSTRNG(20,30) 

20 

1 

NINSTM 

NINSTl 

QST(20) 

STli:N(20) 

KSTRNG(20,30) 

10 

MXFLOW 

NINC 

NINCI 

1 

FLPAR(10,10) 

IFLCOD(IO) 

10 

NINC 

FLPAR(10,10) 

20 

MXTFUN 

NINT 

NINTMI 

ITFCOD(20) 

TPARM(20,20) 

20 

NINT 

TPA1<M<20,20) 

3 

NINT 

NINTMI 

FTF(3) 

TMFRt'O) 

ILF(3) 

Table  4.2.  Problem  Size  Restrictionr 


Restriction 

Limit 

Comments 

Ships/platforms  in  FREQ  SAO 

1 

Program  logic  limit,  requires 
new  coding  to  remove 

Catenary  lines  of  generated 
nodes 

20 

Array  dimensions  in  INPT, 
CATGEN,  INTLIZ  subroutines 

Lines  connected  to  a  node 
where  anchor  holding  power 

CHEK  is  made 

20 

Array  dimensions  and  counter 
check  in  COMCHK  subroutine 

PROB  +  REST  data  sets  in  any 
run 

50 

Array  dimensions  on  title 
cards  read  in  SEADYN  and 

FREINP  subroutines 

Rigid  format  data  sets 

1 

Program  logic  limitation  in 
FREINP  and  SHIPIN  subroutines 

Wave  headings  on  ship  motion 
file 

30 

File  format 

Wave  length  on  ship  motion 
file 

30 

File  format 

Roll  angles  on  ijhip  motion 
file 

8 

File  format 

Wind  velocities  on  ship 
motion  file 

5 

File  format 

Wind  headings  on  ship  load 
file 

20 

File  format 

Current  velocities  on  ship 
load  file 

5 

File  format 

Current  headings  on  ship 
load  file 

20 

File  format 

The  total  size  of  the  CDC  version  of  SEADYN  is  156,000  (octal) 
plus  the  size  required  by  the  ACOM  COMMON.  Selecting  NCOM  to  be  8500 
(decimal),  the  total  program  size  is  176,000  (octal).  A  reduction  in 
total  storage  can  be  achieved  by  segmenting  the  program.  Input  for  one 
approach  using  the  CDC  segmented  loader  is  shown  below. 

TREE  SEAPAY-(INPT,  SAORD,  TRANS,  EIGNS,  FRQSLN,  COMCHK) 

GLOBAL  BUOYS,  CABLE,  CDCAL,  CHKDAT,  COMPNT,  CONTRL,  DSPONI 
GLOBAL  PAYOUT,  RETAPE,  SHIPS,  STRUM,  HEDDAT,  TAPKo,  SIZE,  TYPES 


C.lOBAh 

T BUOYS, 

ICABl.K,  ICNTKL,  IDSBCN, 

GLOBAL 

ISIIIBS, 

1 STRUM,  ITIMKD 

Gi,0BAl. 

liKDCHK . 

LOGIC,  ACOM 

Gl.ORAL 

SHBTAl’, 

SHPLBL,  iFLAG 

FRQSIN 

GLOBAL 

FRQDAT 

KND 

This  approach  to  segmentation  gives  a  maxinuuft  siijc  of  126,000 
(octal)  plus  ACOM.  More  complex  tree  «t nictures  can  be  (b^visefl  to 
shrink  this  size  but  the  complicating  factor  is  th.^it  tin*  FRQSj  N'  biancb 
of  the  tree  calls  the  STATIC  group  of  subroutines  for  wave  Jriit  fovco 
updates  (see  Figure  2.2) 


5.0  DATA  FILES  AND  AUXILIARY  STORAGE 


SEADYN  uses  the  following  FORTRAN  files: 


FILE 

CODE 

Fiu: 

NAME 

MODE 

DESCRIPTION 

01 

bina ry 

Restart  file  for  DEAD  SAO  output;  may  be  used 
for  input. 

02 

binary 

Restart  file  for  LIVE/TSSS  SAO  output;  may  be 
used  for  input. 

03 

binary 

Restart  file  for  DYN  SAO  output;  may  be  used 
for  input. 

04 

binary 

Optional  restart  input  file. 

05 

NIN 

coded 

System  input  file;  read  only  by  FREINP.  Defined 
in  SEADYN  and  carried  in  COMMON /TAPES / , 

06 

NOUT 

coded 

System  output  file;  written  by  various  routines, 
primarily  those  calling  PAGHED  (see  Secti ->0  3.0). 
Defined  in  SEADYN  and  carried  in  COMMON/')  APES/ . 

08 

NSRDC 

binary 

Ship  motion  file;  read  by  SHPMOF  and  SHPRED 
subroutines . 

09 

NTAPE9 

binary 

Scratch  file  for  temporary  storage  of  the  content 
of  COMMON/ACOK/ ,  etc.  in  MODE  and  FREQ  SAOs. 

10 

binary 

Ship  load  file;  written  by  SHIPIN  and  read  by 
SHPIJIS . 

11 

binary 

FRF.Q  steady  state  response  solutions;  written 
by  FRQSLV  and  read  by  FRQSLV. 

1 2 

b  i  n  a  ry 

FREQ  RAO  outputs;  written  by  FRQSLV  and  read  by 
FRQRND  and  FRQRKG. 

13 

NTPCHK 

b  i  ria  ry 

Used  to  pass  the  contents  of  COMMON/ACOM/  from 
FRQSLN  to  CONCHK  for  multiple  wave  heading 
solutions. 

15 

NTAPEl 

b  i  iia  ry 

Deciphered  card  images  from  free-form  input; 
written  by  FREINP  and  read  by  FFLD. 

16 

NTAPF,2 

coded 

Scratch  tile  for  rigid  format  input  data; 
written  by  FREINP  and  read  by  SHIPIN. 

20 

binary 

Mode  shape  (eigenvectors)  output;  written  by 
F.IGNS. 

22 


FILE 

FILE 

CODE 

NAME 

MODE 

DESCRIPTION 

21 

bina ry 

Scratch  file  for  mode  shape 
to  move  the  system  stiffness 
position  in  COMMON/ACOM/  in 
Done  in  EIGNS. 

calculations.  Used 
matrix  to  a  new 
the  needed  format. 

5.1  Restart  File  Structure 


The  SEADYN  program  creates  up  to  three  restart  files  (one  each  for 
the  DEAD,  LIVE,  and  DYN  subanalyses).  Multiple  selections  of  a  subanalysis 
type  simply  extend  the  file  unless  a  rewind  is  signalled  by  the  SAVE 
data  record.  A  counter  (IFILE  array)  is  provided  for  each  of  the  files 
to  keep  track  of  how  many  restart  records  have  been  written.  Each  time 
the  file  is  rewound,  the  counter  for  that  file  is  set  to  zero,  and  a 
label  record  is  written.  The  write  statement  is: 

WRITE (Kl'ILE)  (TITLE(I)  ,1=1  ,NHED)  ,NINA,NPRECZ 
Each  restart  save  operation  uses  the  following  write  statement: 

WRITE(NTAPE)  (A(I ), 1=1 , NINA) , (B(I ) ,1=1 ,NINB) , (C(I) ,1=1 ,NINC) , 

.  (RL(I)  ,I=1,NINRL),(P(I),I=1,NINP0),(T(I),I=1,NINT), 

.  (SH(I) , 1=1 ,NINSHP) , (DP(I) ,I=1,NINDSP) , (STM(I) ,1=1 ,NINSTM) ,NFILE 
.  ,  (I.AABU(I)  ,  1=1  ,NINBI) ,  (lAACAd ) ,  1=1  ,NINCI) 

.  ,  (lAACL(T) ,1=1 .NINCLI) , (lAADS(I) ,I=1,NINDSI) 

.  , (lAAPO(I) ,1  =  1 ,NINPOI) , (lAASHP(l)  ,I=1,N1NSH1) 

.  ,(IAASTM(I),I=1,NINSTI),(IAATIM(I),I=1,NINTMI) 

.  ,  DLD ,  WLD ,  DYN ,  ClIECKR ,  NOVEL ,  NOITER ,  NOFLUD ,  NOLOAD 
.  ,FEEDBK,P0UT,REF1JP,STEPUP 

5.2  ^ Moti on^  File  Structure 

The  Ship  Motion  File  i.s  organized  in  logical  reocrds.  The  specific 
cor, tents  oi  each  record  will  be  described  below.  There  are  seven  distinct 
record  types.  The  first  two  records  contain  data  which  are  independent 
of  wrive  heading  or  wavelength.  Record  types  3  through  7  are  dependent 
on  heading  and  wavelength  and  are  repeated  in  a  nested  loop  fashion. 

Tiie  overall  form  is: 


Record  1 
Record  2 

Looj) 

Loop 

-  Record  3 

on 

on 

Record  4  — 

—  (present  on 

Wave 

Wave 

Record  5 

first  wavelength 

Length 

Heading 

Record  6 

only  for  NCEL 

! 

- Record  7 

formatted  files) 

The  wave  headings  are  assumed  to  be  listed  in  decreasing  order  with  +180 
degrees  being  the  largest  allowed.  The  interpolation  routines  assume 
the  values  given  for  +180  degrees  will  be  used  for  -180  degrees;  therefore, 
data  for  -180  degrees  need  not  be  given. 
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The  wavelengths  are  assumed  to  be  listed  in  decreasing  order  (i.e., 
increasing  frequency  order). 

The  individual  records  of  the  file  are  described  in  terms  of  the 
Fortran  read/write  lists  associated  with  each  record. 

RECORD  1  NA.'ffi  1,  NAME  2,  NAME  3 

Three  Hollerith  variables  providing  identifying  data. 


RECORD  2  (TITO(I),  1=1,  12),  WORD,  WORD  2,  WORD  3,  ELL,  BEAM,  DRAFT,  TVOL, 

TMAS,  TPST,  ZC,  CBV,  NOB,  (FN(I),  1=1,  NOB),  NOH,  (HDGl(I),  1=1,  NOH) , 
NOK,  fBAM(I),  1=1,  NOK),  VNY,  GRAV,  NRV,  (RANG(I),  1=1,  NRV) , 
((GMU(1,J),  J=l,  6),  1=1,6),  ((DC(I,J),  J=l,6),  1=1,6) 


TITO 

WORD 

WORD  2 

WORD  3 

ELL 

BEAM 

DRAFT 

TVOL 

TMAS 

TPST 

ZG 

CBV 

NOB 

FN(I) 

NOH 

HDGI(I) 


NOK 

BAM(I) 

VNY 

GRAV 

NRV 

RANGE(l) 

GMU(I,J) 

DC(I,J) 


=  Hollerith  title  consisting  of  12  6-character 
words 

=  length  unit  label  (6-character  Hollerith) 

=  force  unit  label  (6-character  Hollerith) 

=  moment  «mit  label  (6-character  Hollerith) 

=  Ship ' s  length  (L) 

=  Beam  (L) 

=  Draft  (L)  „ 

=  Ship's  volume;  obtained  from  (ELL/2)'^  .  TVOL 
=  Ship’s  mass  (FtV  ) 

=  Longitudinal  distance  from  c.g.  to  forward  most 
station;  obtained  from  (ELL/2)  .  TPST 
=  Vertical  distance  troro  water  line  to  c.g., 

(+  up)  (L) 

=  Vertical  distance  from  water  line  to  center  of 
buoyancy  (+  up);  obtained  from  ELL  .  CBV 
=  Number  of  speeds  (SEADYN  expects  only  one) 

=  The  Froude  numbers  for  each  speed 
=  Number  of  wave  headings 

=  The  wave  headings  listed  in  decreasing  order 
starting  with  180  degrees  and  proceeding  no 
further  than  -180  degrees  A6 
=  Ntunber  of  wavelengths 

=  Nondimensionai  wavelength  in  decreasing  order, 

\  =  ELL  .  BAM(I)  ,, 

=  Fluid  viscosity  (L''T) 

=  Gravitational  acceleration  (LT  ) 

=  Number  of  roll  angles 

=  The  values  of  roll  anges  (radians)  listed  in 
increasing  order 

=  The  nondimensionai  mass  matrix 
=  The  nondimensionai  hydrostatic  restoring  matrix 


RECORD  3  MM,  HDGI (MM) ,  JJ,  FN(JJ) ,  LL,  BAM(LL) 

MM  =  Heading  number 

JJ  =  Speed  number 

LL  =  Wavelength  niunber 


RECORD  4  ((DA(I,J),  J=l,6),  1=1,6),  ((DB(I,J),  J=l,6),  1=1,6) 

DA(I,J)  =  The  nondimens ional  added  mass  matrix  for  the 
current  combination  of  heading,  speed,  and 
wavelength 

DB{I,J)  =  The  nondimensional  wave  damping  matrix 

RECORD  5  (BOD(I),  B0D(I+3),  BEV(I),  BEV(I+3),  1=1,3) 

BOD,  BEV  =  The  nondimensional  wave  force  coefficient 

RECORD  6  B44S(I),  1=1,  NRV) 

B44S(I)  =  The  nonlinear  roll  dampi,ng  terms  which  are 

added  to  the  damping  matrix  depending  on  the 
size  of  the  roll  angle  (nondimensional) 

RECORD  -/  ((TX(I,J),  1=1,7),  J=l,7),  ((TY(I,J),  1=1,7),  J=l,7), 

((TMI,J),  1=1,7),  J=l,7),  (TP(I),  1=1,7) 

The  complex  nondimensional  drift  force  coefficients 

5 . 3  Ship  Load  File  Format 

The  ship  load  file  is  set  up  to  be  a  library  of  ship  characteristics 
and  loading  tables.  The  file  contains  one  logical  record  for  each  ship 
catalogued.  It  is  written  in  binary  form  with  the  following  FORTRAN 
statement : 

WRITE(IO)  NWIND,  NETHETW,  WNDVEL,  WNPHED,  WNDCOE,  SCALE,  NCRNT, 

NTHETC,  CURVED,  CURHED,  CURCOE,  TDEPTH,  TBLOCK,  TSLT, 

TSAE,  TSAS,  TSWL,  TSB,  TSD,  TSDSP,  TSAP,  SHPCAP,  WLBL, 

CLBL,  LLBL,  VLBL 

Consiiit  Section  8.0  for  descriptions  of  the  variable  and  Section  3.2  for 
array  .sizes. 

5 . 4  Mode  Shape  F i le  Format 

The  output  of  natural  frequencies  and  mode  shapes  can  be  output  on 
file  20.  The  file  contains  one  record  for  each  natural  frequency.  Each 
record  is  written  as; 

WRITE(20)  I,  EIGEN(I),  FREQ,  PERD,  (TGK(J,I),  J=1,NF3) 

where;  I  =  The  mode  number  V 

EIG£N(I)  =  The  circular  frequency  (radians/sec) 

FREQ  =  The  natural  frequency  (Hz) 

PEm  =  The  natural  period 

TGK  =  The  eigenvector  (mode  shape) 

NF3  =  The  number  of  degrees  of  freedom 
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6.0  RESTART  CAPABILITIES 


The  restart  files  mentioned  previously  allow  for  the  recovery  of 
data  for  various  configurations  of  the  modeled  system.  The  methods  for 
using  this  capability  are  explained  in  Section  6.0  of  the  User's  Manual. 

When  restart  is  used  in  a  separate  run  from  the  one  that  created  the 
saved  file,  the  file  may  be  designated  as  any  file  code  recognized  by 
the  system  since  the  file  code  is  given  in  the  input  data.  The  usual 
procedure  is  to  use  codes  01  through  04,  as  explained  in  Section  5.  It 
should  be  recognized  that  normal  execution  may  require  the  re-use  of 
some  files.  By  using  file  code  04  as  restart  input,  the  user  is  guaranteed 
the  file  will  not  be  used  for  output. 

The  restart  procedure  reade  the  title  record  from  the  restart  file 
and  performs  a  comparison  on  the  first  word  if  requested.  It  then  spaces 
down  the  file  to  the  requested  record  number  by  reading  each  record.  The 
last  record  read  then  occupies  the  COMMON  blocks  mentioned  in  Sections  3.1 
and  3.2.  The  restart  activity  is  aborted  if  the  first  word  comparison  fails 
or  if  the  current  size  of  COMMON/ACOH/  is  too  small  to  receive  the  number 
of  words  written  when  the  file  was  created. 


7.0  DESCRIPTIONS  OF  THE  SUBROUTINES 

Each  of  the  subroutines  is  described  in  this  section.  Cross 
references  for  CALLS  and  CALLED  BY  are  also  given.  These  cross 
references  onit  system  library  functions  and  subroutines,  such  as  SQRT, 

COS,  SIN,  etc.  However,  those  routines  calling  EOF,  DATE,  TIME,  and 
ENCODE/DECODE  are  indicated  since  this  may  be  useful  in  conversion  to 
other  nechaines.  These  particular  routines  are  listed  here  for  convenience: 

CALLS  SUBROUTINE 

EOF  FREINP,  FFLD 

DATE  SEADYN  (results  output  by  PACKED} 

TIME  SEADYN  (results  output  by  PACKED) 

ENCODE/DECODE  FREINP,  SELCTS,  SELECT 
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r  NAME 

NO.  OF 
ARCS 

CALLED 

BY 

CALLS 

COMMON 

USED 

DESCRiniON 

ADJMAS 

9 

EQTION 

DYNFBK 

NRAPIT 

MNPRDT 

BUOYS 

CABLE 

CONTRL 

LOGIC 

TAPES 

TIMED 

Adjust  nodal  force  siuonation  for 
tangential  added  mass  (i.e., 
remove  residual  added  mass  from 
tangential  direction) 

AMPHAZ 

5 

FREQREG 

FRQSLV 

SPONLY 

Get  response  amplitude  and  phase 
for  a  set  of  complex  numbers. 

BANDES 

2 

INPT 

Estimate  bandwidth 

BNMOVE 

10 

STATIC 

STPDYN 

NRAPIT 

TVARY 

BUOYS 

CABLE 

CONTRL 

DSPCON 

LOGIC 

PAYOUT 

TIMED 

Impose  boundary  movement  for  MOVE, 
SBUO,  and  fixed  nodes  for  DYN  and 
TSSS. 

BUODYN 

1 

5 

FRQSLV 

CVCTRN 

SPONLY 

TATT 

TRNSHP 

BUOYS 

CONTRL 

FRQDAT 

SHIPS 

TAPES 

Compute  spherical  buoy  dynamic 
coefficients  for  frequency  domain 
equations . 

BUOSTF 

4 

STFGEN 

BUOYS 

Get  stiffness  matrix  for  mooring 
buoys . 

CABTRN 

4 

SHPFIX 

CROSS 

Get  local-to-global  transformation 
matrix  for  line  elements. 

;  lATFRC 

9 

CATSTF 

CATRES 

PCAFX2 

TAPES 

Computer  end  forces  and  stretched 
length  for  bottom-limited  catenary 
element. 

j  CATGEN 

cv.i 
.  1 

17 

INPNOD 

Generate  nodes  along  a  catenary 
curve . 

i  CATRES 

19 

RESIDL 

CATPRC 

CROSS 

MVPRDT 

BUOYS 

CABLE 

CONTRL 

Compute  internal  forces  for  bottom 
limited  catenary  element  for  force 
residual . 

CATSTF 

lA 

STFGEN 

CATFRC 

TATT 

BUOYS 

CABLE 

General  stiffness  matrix  for 
bottom-limited  catenary  element. 

CONTRL 

TAPES 


28 


NAME 

NO.  OF 
ARCS 

CALIJID 

BY 

CALLS 

COMMON 

USED 

DESCRIPTION 

CBLIMP 

12 

FRQSLV 

CABTRN 

CXSLVR 

DRACO 

TATT 

CABLE 

CONTRL 

Calculate  linearized  damping  terns 
for  line  elements  and  assemble  into 
complex  eyaten  matrix. 

CCOMIT 

6 

FRQSLV 

Fora  complex  aiatrix  triple  product 
and  dot  product. 

CD 

1 

STRUMR 

SIMP 

Function  to  conpute  local  drag 
coefficient  ratio  for  strum  cal' 
culation. 

COMBAN 

8 

FRQSLV 

SPONLY 

Complex  simultaneous  equation 
solution  using  syssaetric  banded 
row  storage  format. 

COMCHK 

18 

MAMIPR 

FFLD 

LOUPD 

PAGHED 

RESTART 

BUOYS 

CABLE 

CHKDAT 

COMPNT 

CONTRL 

DSPCON 

HEDCHR 

KEDDAT 

SHIPS 

SIZE 

TAPES 

Check  for  component  adequacy. 

COMPOT 

1 

INPUT 

COMPNT 

Define,  convert,  and  print  con' 
ponent  inventory. 

COMPRP 

13 

UPDATE 

RESIDE 

UPDATE 

DMPFRC 

MPROP 

CABLE 

CONTRL 

PAYOUT 

Computer  element  material  properties, 
strain  and  load  for  multiaaterial 
payout  element. 

CONCNT 

INPT 

BUOYS 

CONTRL 

Count  connections  to  limited  nodes 
and  bodies. 

COMMAS 

FRQSLV 

Get  consistent  mass  matrix. 

CPLOTR 

A1 

MAMIPR 

Dummy  routine  for  eventual  plot 
interface. 

CRCALL 

4 

SHPOEF 

Calculate  ship's  hull  resistance 
coefficient  for  analytical  load 
functions . 

CROSS 

3 

CABTRN 

CATRES 

STRUMR 

UPDATE 

Vector  cross  product. 
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NO.  OF 

CALLED 

COMMON 

ARCS 

BY 

CALLS 

USED 

DESCRIPTION 

CURREN 

2 

LDUPD 

CURUSR 

CABLE 

Define  components  of  current  for 

MANIPR 

TVARY 

CONTRL 

each  node. 

STRUMR 

TIMED 

LOGIC 

CURUSR 

6 

CURREN 

User  subroutine  for  defining 
current  field. 

CVCTRN 

8 

FRQREG 

MLTPLY 

Coordinate  transformation  on  a 

BUODYN 

FRQSLV 

TMLTPLY 

complex  vector. 

CXSLUP 

5 

FRQSLV 

MVPRDT 

Complex  slave  update. 

CXSLVR 

11 

CBLDMP 

MLTPLY 

Store  complex  slave  partition.s 
for  stiffness  or  mass  matrix. 

DEPCOR 

3 

SHPDEF 

Function  to  compute  depth  correc¬ 
tion  for  ship  load  tables, 

DMPFRC 

10 

COMPRP 

CABLE 

Computes  damping  forces  for 

RESIDL 

CONTRL 

material  damping. 

UPDATE 

TIMED 

DRAGCO 

7 

CBLOMP 

User  subroutine  for  computing  drag 

LDUPD 

coefficients . 

DTCALC 

6 

TRANS 

PAGHED 

CABLE 

Estimates  upper  bound  on  time  step 

STPDYN 

CONTRL 

for  DIM  dynamic  .solution. 

SAORD 

HEDDAT 

TAPES 

TIMED 

aYNEBK 

42 

STPDYN 

ADJMAS 

BUOYS 

Transient  solution  of  residual 

NRAPIT 

LOUPD 

CABLE 

feedback  equations  using  Newmark 

MSTRMS 

CONTRL 

Beta  Method. 

RESIDL 

SHIPS 

SHPFIX 

LOGIC 

SLAVED 

TAPES 

STFGEN 

SLVBAN 

TIMED 

EIGNS 

14 

MODES 

STRUMR 

JACOBI 

ACOM 

Set  up  eigenvalue  solution 

EQTION 

41 

STPDYN 

ADJMAS 

BUOYS 

Form  terms  of  dynamic  equation  for 

LDUPD 

CONTRL 

SLl  and  DIM  solutions. 

MSTRMS 

LOGIC 

RESIDL 

TAPES 

SHPFIX 

TIMED 

SI^VLD 

STFMLT 


NAME 

NO,  OF 
ARCS 

CAIJ.B1) 

BY 

CALLS 

COMMON 

USED 

DESCRIPTION 

ERROR 

k 

FREINP 

FFLD 

Error  message  processor  for  free- 
field  reader. 

FFLD 

5 

SEADYN 

INPT 

SAORO 

FRQFRD 

COMCHK 

FRQREG 

FRQRND 

MANIPR 

ERROR 

EOF 

Reads  binary  records  processed 
by  the  free-field  reader. 

FREINP 

5 

SEADYN 

INPT 

SAORD 

FRQFRD 

COMCHK 

FRQREG 

FRQRND 

MANIPR 

ERROR 

EOF 

Reads  binary  records  processed 
by  the  frea-field  reader. 

FREINP 

5 

SEADYN 

EOF 

ERROR 

SELCTS 

SELECT 

SHIFT 

ENCODE 

DECODE 

TAPES 

Fiee- field  reader  —  decodes  free 
field  record.s  and  writes  a  binary 
interpretation  of  the  input  data. 
It  also  copies  rigid  format  data 
over  to  a  scratch  file. 

FKQFRI) 

h 

FRQSLN 

PAGHED 

FFLD 

CONTRL 

DSPCON 

FRQDAT 

HEDDAT 

TAPES 

Read  input  data  for  FREQ  SAO. 

FRQREG 

8 

FRQSLN 

AMPHAIl’; 

CVCTRN 

FFLD 

PAGHED 

CONTRL 

FRQDAT 

HEDDAT 

HEDCHR 

SHIPS 

TAPES 

TYPES 

IFLAG 

Calculate  regular  wave  response 
data  for  FREQ  SAO. 
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N/^ 

FRQRND 


KRQSI.N 


FRQSLV 


INPBOY 

INPBLC 


NO.  OF 

CALLED 

COMMON 

ARGS 

C^AljLS 

_UI^_ 

10 

FRQSLN 

FFLD 

CABLE 

PACKED 

CONTRL 

FRQDAT 

HEDDAT 

HEDCHR 

SHIPS 

TAPES 

TYPES 

IFLAG 

48 

MANIPR 

RESTART 

ACOM 

STATIC 

BUOYS 

UPDATE 

CABLE 

CONTRL 

SHIPS 

TIMED 

TYPES 

IFLAG 

FRQDAT 

HEDDAT 

HEDCHR 

TAPES 

50 

FRQSLN 

AMPHA2 

BUOYS 

BUODYN 

CONTRL 

CBLDMR 

FRQDAT 

CCOMLT 

HEDDAT 

COMBAN 

HEDCHR 

CONMAS 

SHIPS 

CVCTRN 

CXSLUP 

PACKED 

SHPRi:P 

SFONLY 

STFGEN 

TATT 

TENRSP 

TAPES 

ry 

INPl' 

PACKED 

BUOYS 

CONTRL 

HEDDAT 

HEDCHR 

TAPES 

4 

INPT 

PACKED 

BUOYS 

CONTRL 

HEDDAT 

HEDCHR 

TAPES 

DESCRIPTION 

Calculate  randon  response  statist 
tics  for  P'REQ  SAO. 


Controller  for  the  FREQ  SAO. 


Solve  for  frequency  reponses  for 
FREQ  SAO  and  set  up  RAO's. 


Input  body  tables  data. 

Input  body  location  data. 
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NAME 

NO.  OF 
ARCS 

CALLED 

BY 

CAUS 

COMMON 

USED 

DESCRIPTION 

INPGLT 

10 

INPT 

PACKED 

CONTRL 

HEODAT 

HEDCHK 

TAPES 

Input  <;l«aicnt  data  and  generate 
»issing  data. 

INPFLD 

U 

INPT 

PACKED 

CONTRL 

HEDDAT 

HEDCHK 

TAPES 

Input  flow  tables. 

INPFLll 

3 

INPT 

PACKED 

CONTRL 

HEDDAT 

HRDCIIR 

TAPES 

Input  fluid  data 

INPllM 

3 

INPT 

PACKED 

BUOYS 

CONTRL 

HEDDAT 

HEDCIIR 

TAPES 

Input  limit  set  tables. 

INPI.LC 

3 

INPT 

PACKED 

BUOYS 

CONTRL 

HEDDAT 

HEDCHK 

TAPES 

Input  limit  location  data. 

INPMAT 

2 

INPT 

PACKED 

CABLE 

CONTRL 

HEDDAT 

HEDCHK 

LOGIC 

TAPES 

Input  material  tables. 

1  Nl’NOl) 

10 

INPT 

PACKED 

CATGEN 

CONTRL 

HEDDAT 

HEDCHK 

SHIPS 

TAPES 

Input  nodes  and  lines  of  nodes. 

INPSTM 

3 

INPT 

PACKED 

CONTRL 

HEDDAT 

HEDCHK 

TAPES 

STRUM 

Input  strum  string  data. 
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if 


NO.  OE' 

CALLED 

COMMON 

NAME 

ARCS 

BY 

CALLS 

USED 

DESCRIPTION 

INPT 

42 

MANIPR 

BANDRS 

BUOYS 

Input  new  problem  data  and 

COMPDT 

CABLE 

lize  parameters,  calculate 

CONCNT 

CONTRI. 

added  masses,  weight,  etc. 

FF1,D 

PAYOUT 

INPBDY 

SHIPS 

INPBLC 

SIRUM 

INPELT 

TIMED 

INPF1.0 

HEDDAT 

INPFLU 

HF^nCHR 

INPLIM 

TAPES 

1NPLI.C 

TYPES 

INPMAT 

IFLAG 

m 


INl’NOD 

INFSTM 

INPTEN 

INPriiM 

MPROP 

PAGHEP 

SHIPIN 

SHPINP 

INTUZ 


INPTKN  I* 


INPT  PAGHED 


CONIRI, 

UEDDAT 

HKDCHK 

TAPES 


Input  tension  dntn. 


iKPTlM  U 


INPT  PAGHED 


CONTRI. 

HKDOAT 

HKDOHR 

TAPES 

TIMED 


Input  time  limction  tables. 


INTLlZ  A'j 


UPDATE  BUOYS  Initialize  new  piobVem  data,  com* 
MPROP  CABLE  pnte  lea^tths  and  initial  tension 
CONTRI.  data. 

LOGIC 

SHIPS 

TAPES 


1 NTLSA  4 


SAORD 


BUOYS 

CABU; 

CONTRI. 

CHKDAT 

DSPCON 

PAYOUT 

SHIPS 

STRUM 


Initialize  solution  parameters  .and 
set  defaults  of  subanalysis  options. 
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NO.  OF 
ARCS 

CAll.KO 

CAEW 

COMMON 

USED 

DESCRnn'lON 

TIMED 

TYPES 

I  FUG 
RETAPR 

INTLWM 

12 

iNpr 

PAGHED 

BUOYS 

CABLE 

CONTKL 

LOGIC 

SHIPS 

TAPES 

Compute  masses  and  weights  for 
new  problem. 

JACOBI 

13 

EIGNS 

Solves  the  linear  eigenvalue 
problem  using  the  Jocobi  method. 

LDNTRP 

13 

SUPLDS 

Table  look  up  for  ship's  static 
loads . 

LDUPD 

15 

NKAPIT 

STEP 

EQTION 

C^’NFBK 

COMCHK 

TRANS 

CURKF.N 

DRAGOO 

SHPJJIS 

BUOYS 

CABLE 

CHKDAT 

CONTRL 

i>Ayoi;r 

SHIPS 

TIMED 

LOGIC 

TAPES 

Compute  fluid  drag  loads  on  bodies 
lines  and  ships. 

l.iMCHK 

10 

RKSIUI. 

sSTf'GEN 

TVARY 

BUOYS 

CABLE 

CONTRL 

DSPCON 

SHIPS 

TIMED 

TAPES 

LOGIC 

Check  limit,  conditions  and  fix  or 
free  components  as  needed. 

!iAlS 

N/A 

N/A 

SEADYN 

ACOM 

Main  program.  Establishes  the 
size  of  ACOM, 

MAN  I PK 

42 

SEADYN 

COMCHK 

CPLOTR 

CUKREN 

Fl'I.D 

FRQSLN 

INPT 

MODES 

ooTpirr 

PAGHED 

BUOYS 

CABLE 

CHKDAT 

CONTRL 

DSPCON 

PAYoirr 

SHIPS 

RETAPF. 

STRUM 

Manipulative  routine  to  control 
the  execution  of  the  subana lysis 
modules. 

35 


NAME 


CALLS 


COMMON 

USED  DESCRIFflON 


NAME 

NO.  OF 
ARCS 

CALLED 

BY 

CALLS 

COMMON 

USED 

DESCRIFflON 

SAORD 

STATIC 

STRUMR 

TRANS 

UPDATE 

LOGIC 

HEDDAT 

HEDCHR 

TAPES 

TIMED 

SIZE 

TYPES 

IFLAG 

MASMOD 

13 

PAYUPD 

BUOYS 

CABLE 

CONTRL 

PAYOUT 

SHIPS 

TAPES 

Make  mass  and  weight  adjustments 
tor  payout  elements. 

MLTPLY 

9 

CVCTRN 

CXSLVR 

SHPFIX 

SHPSTF 

SLAVER 

STFGEN 

STRUMR 

TATT 

Matrix  multiply. 

MODES 

90 

MAN I PR 

EIGNS 

STFGEN 

ACOM 

CONTRL 

DSPCON 

Set  up  and  execute  MODE  SAO- 
eigenvalue  solution. 

KPHOP 

10 

COMPRP 

INPMAT 

INPT 

RESIDE 

TENRSP 

UPDATE 

Material  property  table  lookup. 

MSTKMS 

0 

DYNFBK 

EQTION 

NRAPIT 

Dummy  routine  intended  for  event 
ual  use  with  transient  dynamics 
with  slave  nodes. 

MVPRDT 

5 

ADJMAS 

CATRES 

CXSLUP 

SHPLDS 

SHPSTF 

SLAVUP 

STFGEN 

Premultiply  a  vector  by  a  matrix 
and  get  the  dot  product  (square) 
of  the  result. 

NO.  OF 

CALLED 

COMMON 

NAME 

ARGS 

BY 

CALLS 

USED 

DESCRIPTION 

KRAPIT 

41 

STATIC 

ADJMAS 

BUOYS 

Modified  Newton-Raphson  solution 

TRANS 

BNMOVE 

CABLE 

routine  for  statics  and  dynaenics. 

DYNFBK 

CONTRL 

COUPD 

DSPCON 

MSTRMS 

SHIPS 

OUTPUT 

TIMED 

PAGHED 

LOGIC 

RESIDE 

HEDDAT 

SERCHl 

HEDCHR 

SHPFIX 

TAPES 

SLAVED 

SLAVUP 

SLVBAN 

STFGEN 

STFMLT 

TVARY 

OUTPUT 

8 

SEAOYN 

PAGHED 

CONTRL 

Output  routine  for  node  positions 

MANIPR 

RESTART 

PAYOUT 

and  velocities  and  element  tensions 

NRAPIT 

TIMED 

STEP 

HEDDAT 

STATIC 

HEDCHR 

TRANS 

TAPES 

LOGIC 

SIZE 

SHIPS 


PAGHED 


COMCHK 

DTCALC 

FRQFRD 

FRQKEG 

F.KQSLN 

FRQSLV 

IHP13DY 

INPBLC 

INPELT 

INPFLO 

INPFLU 

INPIJM 

INPLLC 

INPMAT 

INPNOD 

INPSTM 

INPT 

INPTEN 

INPTIM 

MANIPR 

NRAPIT 

OUTPUT 

PAYIJPD 


HEDDAT  Write  page  heading  for  output. 

HEDCHK 

SIZE 

TAPES 
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NAME 

NO.  OF 
ARCS 

CALLED 

BY 

CALLS 

COMMON 

USED 

DESCRIPTION 

SAORD 

SEADYN 

SHPINP 

SHPMOF 

STPDYN 

PAYUPD 

21 

STATIC 

STPDYN 

MASMOD 

PAGHED 

TVARY 

CABLE 

CONTRL 

PAYOUT 

TIMED 

TAPES 

HEDDAT 

HEDCHR 

LOGIC 

Compute  pay-out/reel-in  data  for 
changiag  lengths  and  mitosis. 

PCAFX2 

14 

CATFRC 

Compute  end  forces  for  a  catenary 
element  with  given  end  positions. 

PSIN 

1 

STRUMR 

SIMP 

Compute  mode  shape  amplitude  for 
strum  drag  amplification. 

RAND IN 

3 

TVARY 

Dum&y  routine  for  NCEL  random 
signal  input. 

RESIDE 

40 

STEP 

NRAPIT 

EQTION 

DYNFBK 

COMPRP 

DMPFRC 

LINCHK 

MPROP 

SHPSTF 

BUOYS 

CABLE 

CONTRL 

DSPCON 

PAYOUT 

Compute  internal  nodal  point  forces 
(element  forces)  for  the  force 
residual  vector.  Also  updates 
lengths,  tension,  and  material  pro¬ 
perties  for  iterative  solutions. 

SHIPS 

STRUM 

TIMED 

TAPES 

LOGIC 


fiK START  5 


SEADYN 

OUTPUT 

FRQSLN 

COMCHK 


ACOM 

BUOYS 

CABLE 

CONTRL 

DSPCON 

PAYOUT 

TIMED 

SHIPS 

STRUM 

LOGIC 

SIZE 

TAPES 

HEDDAT 

MEDCHK 


Read/write  restart  data  records. 
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NAME 

SAOCUR 

SAOFOF 

SAOIMP 

SAOINI 

SAOLOD 

SAOLVA 

SAOMOV 

SAOPAY 

SAORD 


NO.  OF 
ARCS 

CALLED 

BY 

CALLS 

COMMON 

USED 

3 

SAORD 

HEDDAT 

LOGIC 

TAPES 

5 

SAORD 

PACKED 

CONTRL 

HEDDAT 

TAPES 

0 

SAORD 

BUOYS 

HEDDAT 

TAPES 

6 

SAORD 

HEDDAT 

TAPES 

8 

SAORD 

HEDDAT 

LOGIC 

TAPES 

4 

SAORD 

HEDDAT 

TAPES 

6 

SAORD 

DSPCON 

HEDDAT 

TAPES 

TIMED 

TYPES 

3 

SAORD 

HEDDAT 

LOGIC 

PAYOUT 

TAPES 

DESCRIPTION 

Process  CURR  input  record. 

Process  FIX/FREE  input  records. 

Process  IMPA  input  record. 

Process  INIT  input  records. 
Process  LOAD  input  records. 

Process  LVAR  input  record. 
Process  MOVE  input  records. 

Process  PAYO  input  records. 


9 


MANIPR 


DTCALL 

FFIJD 

INTLSA 

PACKED 

SAOCUR 

SAOFOF 

SAOIMP 

SAOINI 

SAOLOD 

SAOLVA 

SAOMOV 

SAOPAY 

SAOSAV 

SAOSBN 

SAOSOL 

SAOSTP 

SAOSOW 

SAOTIM 


BUOYS 

CABLE 

CONTRL 

CHKDAT 

DSPCON 

HEDDAT 

HEDCHR 

PAYOUT 

SHIPS 

STRUM 

TAPES 

TIMED 

TYPES 

IFLAG 

SIZE 

KETAPK 

LOGIC 


Controlling  routine  of  input  and 
initializing  subanalysis  option 
(SAO)  data. 
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NO.  OF 
ARCS 

CALLED 

BY 

CALLS 

COMMON 

USED 

DESCRIPTION 

SAOSAV 

4 

SAORD 

HEDDAT 

TAPES 

TYPES 

Process  SAVE  input  record. 

SAOSBU 

5 

SAORD 

BUOYS 

HEDDAT 

TAPES 

Process  SBUO  input  record. 

SAOSOL 

3 

SAORD 

CONTRL 

HEDDAT 

TAPES 

TIMED 

Process  SOLU  input  record. 

SAOSOW 

1 

SAORD 

HEDDAT 

SHIPS 

TAPES 

Process  SURF/WIND  input  records. 

SAOSTP 

7 

S.AORD 

HEDDAT 

TAPES 

TYPES 

Process  STEP  input  record. 

SAOTIM 

4 

SAORD 

HEDDAT 

TAPES 

TYPES 

Process  TIME  input  record. 

SEADYN 

1 

MAIN 

DATE 

FFLD 

ACOM 

BUOYS 

Primary  controlling  routine.  Pro¬ 
cesses  title  cards  and  PROB/REST 

FREINP  CABLE  data  records.  Call  for  restart 


MANIPR  CDCAL  read  (if  needed)  and  establishes 
PAGHEO  CHKDAT  the  variable  dimension  position 
RESTART  COMPNT  of  main  data  items  in  ACOM. 
SEAMES  CONTRL 
TIME  DSPCON 
LOGIC 
PAYOUT 
RETAPE 
SHIPS 
STRUM 
TAPES 
TIMED 
TYPES 
IFLAG 
HEDDAT 
HEDCHR 
SIZE 

i 

SEAMES  0  SEADYN  HEDDAT  Prints  user  information  messages 

HEDCHR  at  the  beginning  of  output. 

TAPES 
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NAME 

NO.  OF 
ARCS 

CALLED 

BY 

CALLS 

COHMON 

USED 

DESCRIPTION 

SELCTS 

1 

FREINP 

ENCODE 

Interprets  the  solution  option 
input  codes  and  trsnslates  to 
numeric  value. 

SELECT 

2 

FREINP 

ENCODE 

TYPES 

IFLAG 

Interpretes  the  keywords  from 
free- field  input  and  translates 
to  nuMtic  values. 

SERCHl 

7 

NRAPIT 

Performs  linear  interpolation  for 

ID  search  accelerator  for  MNR 
solution. 

SHI PIN 

0 

INPT 

CONTRL 

SHIPS 

TAPES 

SHPTAP 

SHPLBL 

Reads  the  static  ship  load  file 
data  placed  on  file  code  NTAPE2 
as  Rigid  format  data  by  FREINP 
and  translates  it  to  the  ship  load 
file  (binary  form). 

SHPDEF 

3 

SHPINP 

CRCALC 

DEPCOR 

HEDDAT 

HEDCHR 

SHPLBL 

SHPTAP 

TAPES 

Process  SHIP  input  data. 

SHPFIX 

4 

DYNl'BK 

EQTION 

NRAPIT 

STEP 

CABTBN 

MLTPLY 

TMLPLY 

BUOYS 

SHIPS 

Adjust  residual  for  ship  and  moor'* 
ing  buoy  constraints.  Unless  SHIPK 
data  are  input,  ships  are  fixed  in 
HEAVE,  PITCH,  and  ROLL  for  static 
analyses.  Unless  more  than  two 
line  attachments  are  made  to  moor¬ 
ing  buoys,  their  ROLL  response 
will  be  fixed. 

SHPINP 

3 

INPT 

PACKED 

SHPDEF 

CONTRL 

HEDDAT 

HEDCHR 

SHIPS 

TAPES 

Set  up  for  reading  ship  input 
data. 

SHPLDS 

2 

IDUPD 

LDNTRP 

HVPRDT 

TRNSHP 

CONTRL 

SHIPS 

SHPLBL 

SHPTAP 

Compute  ship  static  loads . 

SHPMOF 

1 

FRQSLN 

PACKED 

CONTRL 

FRQDAT 

Initialize  data  from  ship  motion 
file. 

HEODAT 

HEDCHR 

SHIPS 

TAPES 


NO.  OF 

CALLED 

COMMON 

-.il. 

ARCS 

BY 

CALLS 

USED 

DESCRIPTION 

-  ■ 

SHPRED 

24 

FRQSLV 

Read  data  and  interpolate  from 

ship  motion  file. 

SHPSTF 

9 

STFGEN 

MLTPLY 

SHIPS 

Compute  and  assemble  ship  stiffness 

--3 

KESIDL 

MVPRDT 

matrix  or  force  residual  contribu¬ 

- 

tion. 

'fl 

SIMP 

6 

STRUMR 

F 

Simpson  rule  quadrature.  External 

1 

(PSIN) 

function  F  provided  in  calling 

(CD) 

sequence . 

IMPID 

18 

STFGEN 

MVPRDT 

TAPES 

Compute  stiffness  matrix  for  ID 

'jii 

TATT 

simplex  element  (cable/truss)  and 

, 

assemble  in  global  stiffness 

\r 

Jj 

:7 

matrix. 

■  ■’? 

SLAVER 

10 

STFGEN 

MLTPLY 

Store  stiffness  terms  for  master 

node  rotation. 

••if 

1 

SLAVLD 

7 

DYNFBK 

Transfer  loads  at  slave  nodes  to 

EQTION 

master  nodes. 

NRAPIT 

in 

j 

STEP 

-}ir 

TRANS 

'i 

SLAVUP 

8 

NRAPIT 

MVPRDT 

Recover/ update  displacements  and 

.$ 

STEP 

TRANSHP 

positions  for  slave  nodes. 

UPDATE 

SLVBAN 

8 

DYNF'BK 

Decompose  and  solve  system  of 

'  : 

NRAPIT 

simultaneous  linear  algebraic 

-b 

STEP 

equations.  Assumes  compacted 

i 

symmetiic  banded  row  format. 

1 

SrONLY 

li 

BUODYN 

AMPHAZ 

Solve  for  unrestrained  ship 

FRQSLV 

COMBAN 

response . 

i 

.  1 

STATIC 

46 

MANIPR 

BNMOVE 

CONTRI. 

Control  the  static  solutions: 

FRQSLN 

NRAPIT 

DSPCON 

DEAD,  LIVE,  TSSS. 

■  i 

-  ^  j 

OUTPUT 

HEDDAT 

' 

PAYUPD 

HEDCHR 

STEP 

LOGIC 

1 

STRUMR 

PAYOUT 

:j 

TVARY 

RETAPE 

UPDATE 

SHIPS 

"1 

STRUM 

SIZE 

TAPES 

TIMED 
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NO.  OF 

CALLED 

COMMON 

Nj^ 

ARCS 

BY 

CALLS 

USED 

STEP 

41 

STATIC 

LDUPD 

BUOYS 

OUTPUT 

CABLE 

RESIOL 

CONTRL 

SHPFIX 

DSPCON 

SLAVED 

HEDCHR 

SLAVUP 

HEDDAT 

SLVBAN 

LOGIC 

STEPAD 

SHIPS 

STEPDP 

STEPOV 

TAPES 

STEPAD 

22 

STEP 

SLAVUP 

BUOYS 

CONTRL 

DSPCON 

HEDCHR 

HEDDAT 

TAPES 

STEPDP 

11 

STEP 

BUOYS 

CABLE 

TAPES 

CONTRL 

DSPCON 

TIMED 

LOGIC 

SHIPS 

STEPOV 

12 

STEP 

CONTRL 

DSPCON 

TAPES 

STFGEN 

40 

DYNJBK 

BUOSTF 

BUOYS 

EQTION 

CATSTF 

CABLE 

MODES 

LIMSHK 

CONTRL 

NRAPIT 

MLTPLY 

DSPCON 

STEP 

MVPRDT 

PAYOUT 

STPDYN 

SHPSTF 

SHIPS 

FRQSLV 

SIMPID 

STRUM 

STRUMR 

SLAVER 

TAPES 

TIMED 

STFMLT 

7 

EQTION 

NRAPIT 

STPDYN  41 

TRANS 

BNMOVE 

CONTRL 

DTCALC 

HEDCHR 

DYNFBK 

HEDDAT 

EQTION 

LOGIC 

PACKED 

PAYOUT 

DESCRIPTION 

Static  solution  using  SLI,  RFB, 
VRS,  VRR  methods. 


Makes  solution  parameter  adjustments 
for  the  VRS,  VRR  methods. 


Compute  residual  ADTM  and  damping 
for  VRR  solution. 


Adjust  for  limit  overshoot  for 
STEP  subroutine. 


Generate  global  stiffness  matrix. 


Matrix/vector  multiply  using 
compact  banded  stiffness  matrix 
format . 

Solution  routine  for  DYN  SAO. 
Does  SLI »  HNR,  RFB,  and  D1 
solutions . 
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NO.  OF 

CALLED 

COMMON 

NAME 

ARCS 

BY 

CALLS 

USED 

DESCRIPTION 

PAYUPD 

TAPES 

STFGEN 

TVARY 

TIMED 

STRUHR 

40 

MANIPR 

CD 

BUOYS 

EstiaateB  drag  amplification  due 

STATIC 

CROSS 

CABLE 

to  strunning. 

TRANS 

CURREN 

CONTRL 

EIGNS 

LOGIC 

MLTPLY 

SIZE 

PSIN 

STRUM 

SIMP 

STFGEN 

TATT 

TIMED 

TATT 

8 

CATSTF 

SIMPID 

CBLDMP 

STRUMR 

FRQSLV 

MLTPLY 

Coordinate  transformation  of  matrix 
A  with  matrix  T  using  the  form 

T  A  T^. 

TEN8SP 

11 

FRQSLV 

MPROP 

CABLE 

Compute  components  of  dynamic 
tension  for  FREQ  SAG. 

TFNUSR 

4 

TVARY 

User  written  subroutine  for  time 
functions . 

TMLPLY 

9 

CVCTRN 

SHPFIX 

Pre-multiply  a  mat.|pix  by  a  matrix 
transposed.  C  =  A*  B. 

TRANS 

44 

MANIPR 

LDUPD 

BUOYS 

Controller  for  transient  dynamics 

NRAPIT 

CABLE 

(DYN) 

OUTPUT 

CONTRL 

SIJVVLD 

DSPCON 

STPDYN 

HEDCHR 

STRUMR 

HEDDAT 

TVARY 

LOGIC 

UPDATE 

PAYOUT 

RETAPE 

SKIPS 

SIZE 

STRUM 

TAPES 

TIMED 


I’RANSHP  4 


SHPLDS 

SLAVUP 

UPDATE 


Coapute  local-to'global  transforma¬ 
tion  for  a  given  heading. 


NO.  OF 

CALLF.D 

COMMON 

NAME 

ARCS 

BY 

^LLS 

USED 

DESCRIPTION 

TVARY 

3 

BNMOVE 

RANDIN 

CONTRL 

Compute  time  variatioa  functions. 

CURREN 

LIMCHK 

NRAPIT 

PAYUPD 

STATIC 

STPDYN 

TRANS 

UPDATE 

TFNUSR 

TIMED 

UPDATE 

41 

HANIPR 

COMPRP 

BUOYS 

Update  displacements,  positions, 

INPLIZ 

CROSS 

CABLE 

tensions,  strains,  and  lengths. 

FRQSLN 

DMPFRC 

CONTRL 

Get  new  coordinate  transformation 

STATIC 

MRPROP 

DSPCON 

data  as  needed. 

TRANS 

SUVUP 

LOGIC 

TRNSHP  PAYOUT 
TVARY  SHIPS 
TAPES 
TIMED 
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8.0  DESCRIPTIONS  OK  MAJOR  VARIABLES 


This  section  lists  the  major  variables  in  alphabetical  order  and 
gives  a  brief  description  of  each.  Entries  in  the  ARRAY  SIZE  indicate  the 
dimension  of  the  array.  No  entry  means  a  scalar  variable.  A  variable  name 
(e.g.,  NE)  means  an  array  which  is  mapped  to  the  variable  dimension  working 
space  of  COHMON/ACOM/  through  the  HANIPR  subroutine  calling  sequence.  The 
length  of  such  arrays  are  NE  (number  of  elements)  or  NN3  (three  times  the 
number  of  nodes).  The  entries  under  the  type  column  are; 

R  —  for  real  variable  (floating  point) 

I  —  for  integer  variable 

C  -•*  for  complex  variable 

1.  --  for  logical  variable 

H  --  for  Hollerith  word  for  title  or  label  output 

The  column  named  INPUT  KEYWORD  indicates  the  input  record  which  provides 
the  initial  definition  of  the  variable.  RIGID  refers  to  the  rigid  format 
input  for  the  ship  load  file. 
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VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYj^ 

DESCRIPTION 

A 

NE 

ACOM 

R 

Cable  element  diameters. 

AABUOY 

100 

BUOVS 

R 

Additional  array  for  future 
expans  Ion . 

AACAB 

100 

CABLE 

R 

Additional  array  for  future 
expansion. 

AACON 

100 

CONTRL 

R 

Additional  array  for  future 
expansion. 

AADSP 

100 

DSPCON 

K 

Additional  array  lor  future 
expansion. 

AAPAY 

100 

PAYOUT 

K 

Additional  array  for  future 
expansion. 

AASliP 

100 

SHIPS 

R 

Additional  array  for  future 
expansion. 

AAS'l'M 

100 

STRUM 

H 

Additional  array  for  future 
expansion . 

AATIM 

100 

TIMED 

R 

Additional  array  for  fut  ii't 
expansion. 

AHUOY 

TYPES 

I 

Numeric  code  for  BUOY  keyword 
=  S6 

At.CCUK 

SHIPS 

R 

Acciunul aloil  surface  current 
load  factor  for  DKAD/LIVE. 

A''Ci'Ai’ 

CONTRL 

R 

Accumulated  load  factor  for 
DEAD/1.1  VE. 

ACCWND 

SHIPS 

R 

Accmnulated  wind  load  factor 
for  DKAD/LIVE. 

A(.;ONE 

TYPES 

1 

Nua^eric  code  for  CONE,  keyword 
=  57. 

ADM 

so 

BUOYS 

K 

Buoyant  forces  for  BODY 
table. 

INPUT 

KEYWORD 


BODY 


AUKAI) 


TYPES 


1 


Num«*ric  co'le  for  HEAD  keyword 

=  il. 


VARIABLV. 

NAME 

AKRAY 

SJl  lK  . 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

mww) 

AIMPA 

TYPES 

1 

Niunoric  codr  for  IMPA  kryvord 
=  Tl. 

AIN/f 

TYPES 

1 

Numeric  codr  for  INIT  kryword 
=  32. 

.KVE 

TYPES 

1 

Numeric  code  for  INVE  keyword 
-  18. 

ALIMI 

TYPES 

1 

Niuneric  code  for  LI  MI  keyword 
=  19. 

ALINK, 

TYPES 

I 

Numeric  code  for  LINE  keyword 
=  20. 

ALIVE 

TYPES 

I 

Numeric  code  for  LIVE  keyword 
=  2. 

ALi-OC 

TYPES 

1 

Numeric  code  for  LLOC  keyword 
=  21 

AlOAl) 

TYPES 

1 

Niuneric  code  for  LOAD  keyword 
-  34. 

AI.PNEW 

TIMED 

R 

Generalized  Newmark  iutegra- 
liou  paramelrr. 

TIME 

ALVAK 

TYPES 

I 

Numeric  code  for  I.VAR  keyword 
=  Tf). 

AMATK 

TYPES 

I 

Numeric  code  for  MATE  keyword 
--  22. 

AAtAX 

CDCAI, 

H 

Strum  string  response 
p.iramcter . 

5 

PAYOUT 

K 

Mitosis  lengths  for  payout 
elements . 

PAYO 

AMODF. 

TYPES 

1 

Niuneric  code  for  MODE  keyword 
= 

amove: 

TYPES 

I 

Numeric  code  for  HOVF;  keyword 
=  .36. 

AMPMN 

KRQDAT 

R 

Cut  oil  wave  amplitude  for 

FREQ  SAO. 
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VARIABl.K 

NAME 

ARRAY 

S 1 ZE 

COMMON 

NAME 

TYPE 

DKSCKl P riON 

INPUT 

KEYWORD 

AMSOI. 

TYPES 

1 

Numerit'  coHi‘  for  MSOL  k«?ywoi'<I 
=  3/ . 

ANCII 

TYPES 

1 

Numeric  code  for  ANCH  keyword 
=  55. 

ANCTAB 

5,16,6 

COMPNT 

K 

Inventory  anchor  properly 

table . 

ANEW 

TYPES 

1 

Numeric  code  for  NEW  keyword 
=  9. 

ANODE 

TYPES 

1 

Numeric  code  for  NODE  keyword 
=  23. 

APROP 

3 

SHIPS 

R 

Propeller  projected  area. 

SHIP 

ASH  IP 

TYPES 

I 

Numeric  code  lor  SHIP  keyword 
=  24. 

ASTEP 

TYPES 

I 

Niuneric  code  for  STEP  keyword 
-  43. 

ASTRUM 

TYPES 

1 

Numeric  code  lor  STRU  keyworil 
=  25. 

ATI  ME 

TYPES 

I 

Numeric  code  lor  TIME  keyword 
=  45. 

ATYPE 

3,t> 

COMPNT' 

K 

luvcutory  aiuTior  type  table. 

AWIND 

TYPES 

I 

Numeric  code  for  WIND  keyword 
~  46. 

b/xM 

30 

FRQDAT 

R 

l.i.st  wavelength  from  ship 
motion  file. 

KAMC 

50 

BUOYS 

R 

Added  mass  coefficients  for 
body  table. 

BODY 

BKTNEU 

TIMED 

R 

Integration  parameter,  p. 

TIME 

HERN 

50 

BUOYS 

R 

(.ength  for  cylindrical  buoys, 
body  t.able. 

BODY 

BLOC 

TYPES 

1 

Numeric  code  lor  BI.OC  keyword 
=  13. 

mi 
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VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

BLOCK 

3 

SHIPS 

R 

Block  coefficient  for  each 
ship. 

BMOM 

50 

BUOYS 

R 

Mass  moment  of  inertia  for 
body  table. 

BODY 

BODY 

TYPES 

I 

Numeric  code  for  BODY  keyword 
=  14. 

BOMAS 

50 

BUOYS 

R 

Body  mass  for  body  table. 

BOVOL 

50 

BUOYS 

R 

Body  volume  for  body  table. 

BSCD 

50 

BUOYS 

R 

Surface  current  drag  for 
body  table. 

BODY 

BTYPE 

3,2 

COMPNT 

R 

Inventory  buoy  type  table. 

BUOKP 

3,50 

BUOYS 

R 

Body  location  parameters  for 
checking  static  response 
instability. 

BUOTAB 

7,6,2 

COMPNT 

R 

Inventory  buoy  property 
table. 

BWND 

50 

BUOYS 

R 

Wind  drag  coefficients  for 
surface  buoy. 

BODY 

B44S 

8 

FRQDAT, 

R 

Roll  damping  terms  for  ship 
motion  file. 

CABMAS 

10 

CABLE 

R 

Mass  per  unit  length  for 
cable  materials. 

CAI) 

SHIPS 

R 

Global  current  heading  in 
degrees . 

SURF 

CAMC 

10 

CABLE 

R 

Line  material  added  mass 
coefficient . 

MATE 

CBUO 

BUOYS 

R 

Body  drag  coefficient. 

CEPS 

STRUM 

R 

Strum  update  parameter. 

SOLU 

CHAIN 

31,3 

COMPNT 

R 

Inventory  chain  property 
table . 

CHECKR 

LOGIC 

L 

Component  check  flag;  True 

means  CHEK  SAO  active. 
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VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

CHEK 

TYPES 

I 

Numeric  code  for  CHECK  key¬ 
word  =  7 . 

CLBL 

SHPLBL 

H 

Current  load  label  for  ship¬ 
load  file. 

RIGID 

CORDLM 

50 

BUOYS 

R 

Limit  set  coordinate  for 
limit  table. 

LIMI 

CPROP 

3 

SHIPS 

R 

Propeller  resistance 
coefficient. 

SHIP 

CR 

3 

SHIPS 

R 

Longitudinal  resistance 
coefficient . 

SHIP 

CRST 

30 

STRUM 

R 

Not  used  at  present. 

cs 

3 

SHIPS 

R 

Hull  wetted  surface 
coefficient . 

SHIP 

CURCOE 

20,3,5 

SHIPS 

R 

Current  load  coefficient 
for  ship  load  tables . 

RIGID 

CURBED 

20 

SHIPS 

R 

Current  headings  for  ship 
load  tables. 

RIGID 

CURLEN 

5 

PAYOUT 

R 

Current  length  of  payout 
element . 

CURMUL 

CABLE 

R 

Flow  field  scale  factor. 

CURR 

viURNT 

SHIPS 

R 

Magnitude  of  surface 
current . 

SURF 

CURR 

TYPES 

I 

Numeric  code  fox  CURR  key¬ 
word  =  28. 

CYLRMS 

50 

BUOYS 

R 

Added  mass  for  cylindrical 
buoy  in  body  table. 

DA 

6,6 

FRQDAT 

R 

Ship's  added  mass  matrix  from 
ship  motion  file. 

DALPHA 

TIMED 

R 

Proportional  damping  multi¬ 
plier  of  mass  matrix  or 
alpha  integration  parameter 
for  VRS/VRR  .''.olutions . 

SOLU 

/ 
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VARIABLk'. 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

DB 

6,6 

F8QDAT 

R 

Ship’s  damping  matrix  from 
ship  motion  file. 

DBETA 

TIMED 

R 

Proportional  damping  multi¬ 
plier  of  stiffness  matrix 
or  initial  pseudo-time  step 
for  VRS/VRR  solutions. 

SOLD 

DBU 

50 

BUOYS 

Body  diameters  for  body 
table. 

BODY 

DC 

6 ,6 

FRQDAT 

R 

Ship's  restoring  force 
matrix  from  ship  motion 
file. 

DEAD 

TYPES 

1 

Numeric  code  for  DEAD  key¬ 
word  =  1 . 

DELFAC 

CONTRL 

R 

Load  factor  increment  (step 
size)  for  DEAD/LIVE. 

DELTMP 

CONTRL 

R 

Temporary  storage  for  DELFAC 
for  VRS/VRR  solutions. 

DEPiAD 

CONTRL 

K 

Degree/radian  conversion 
factor. 

DERR 

CONTRL 

R 

Displacement  error  tolerance 
for  iterative  solutions. 

SOLU 

DIArt 

10 

CABLE 

Line  material  diameters. 

MATE 

DISPC 

:30 

DSPCON 

R 

Static  imposed  displace¬ 
ment  amplitudes  in  C. 

MOVE 

DISPP 

30 

DSPCON 

R 

Rcf'-ence  imposed  displace¬ 
ment  amplitudes. 

DLD 

LOGIC 

L 

Dead  load  flag;  True  means 

DEAD  SAO. 

DLREF 

PAYOUT 

R 

Unstretched  lengths  for  pay¬ 
out  elements. 

DMAXAB 

TIMED 

R 

Magnitude  of  maximum  accele¬ 
ration  increment  in  DYN  SAO. 

DMAXP 

TIMED 

K 

Previous  maximum  acceleration 

in  DYH  SAO. 
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VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME  TYPE 

DESCRIPTION 

INPUT 

DMD 

KEYWORD 

CONTRL  R 

Numerical  damping  coeffi¬ 
cient. 

SOLU 

DOMG 

FRQDAT  R 

Wave  frequency  increment 

SPEC 

DONE 

TYPES  I 

Numeric  code  for  DONE  key¬ 
word  =  54. 

DPTMF 

5 

PAYOUT  R 

Time  function  increment  for 
payout  point. 

DRAD 

SNIPS  R 

Degrees  per  radian. 

DRFTFR 

3 

FRQDAT  R 

Wave-induced  drift  force 
components . 

DRGAMP 

NE 

ACOM  R 

Drag  amplification  factor. 

DS 

NE 

ACOM  R 

Current  lengths  for  ele- 
meutu . 

DSO 

NE 

ACOM  R 

Initial  lengths  for  ele¬ 
ments. 

ELEM 

DSR 

NE 

ACOM  R 

Reference  lengths  for  ele¬ 
ments  . 

DT 

TIMED  R 

Time  step. 

TIME 

OTH 

TIMED  R 

Upper  bound  on  time  step. 

DTL 

TIMED  R 

Lower  bound  on  time  step. 

DTLL 

TIMED  R 

Lower  bound  on  time  step  in 
DYN  SAO. 

DTMAX 

TIMED  R 

Upper  bound  on  time  step  in 
DYN  SAO. 

DTRSRT 

TIMED  R 

Restart  save  time  interval. 

SAVE 

DTU 

TIMED  R 

Reference  update  time. 

TIME 

DU 

NE 

ACOM  R 

DYN 

LOGIC  L 

Transient  dynamic  flag;  True 

■leans  DYN  SAO. 


VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

E 

2,10 

CABLE 

R 

Slopes  of  line  segments  in 
material  tables. 

. 

4 

ELEM 

TYPES 

I 

Numeric  code  for  ELEM  key¬ 
word  =  15. 

ELL 

FRQDAT 

R 

Ship's  length  for  nondimen- 
sionalized  ship  motion  file. 

ERR 

ERR 

R 

Not  used. 

. 

1 

ES 

NE 

ACOM 

R 

Element  secant  EA  in  ^C. 

ESP 

2,5 

PAYOUT 

R 

Multimaterial  payout  secant 
modulii. 

ET 

NE 

ACOM 

R 

Element  tangent  EA  in  ^C. 

s 

Em 

TYPES 

I 

Numeric  code  for  EXTE  key¬ 
word  =  49. 

1 

EXTRAP 

CONTRL 

R 

MNR  extrapolation  factor. 

F 

NN3 

ACOM 

R 

Nodal  loads  in  ^C. 

FACH 

SHIPS 

R 

Depth  correction  factor. 

FACLEN 

FRQDAT 

R 

Length  conversion  factor 
for  ship  motion  file. 

EXTE 

FBUSk 

50 

BUOYS 

R 

Surface  buoy  time  function 
value  for  body  locations. 

FDEPTH 

5 

CABLE 

K 

F'luid  interface  depths. 

FLUI 

1 

-L” 

FDIVY 

CONTRL 

R 

Inventory  conversion  factor 
for  line  diameters. 

INVE 

FEEUBK 

LOGIC 

L 

F'eedback  flag;  True  means 
feedback  form  used,  RFB  or 
VRR. 

FELT 

3,2 

CHKDAT 

R 

Fluid  loads  on  a  single 
element  for  CHEK  SAO. 

FG 

NN3 

ACOM 

K 

Nodal  point  gravity  loads. 

1 

FGAM 

5 

CABLE 

R 

Fluid  table  specific 

FLUI 

1 

weights . 
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VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT  ^ 

KEYWORD  ^ 

FINC 

CONTRL 

R 

Alters  the  residual  added 
mass  correction  for  TL  or 

UL  solutions. 

FINIS 

TYPES 

I 

Numeric  code  for  END  key¬ 
word  =  10. 

-li 

FIX 

TYPES 

I 

Numeric  code  for  FIX  key¬ 
word  =  29. 

FLNVY 

CONTRL 

R 

Inventory  conversion  factor 
for  buoy  diameters  and 
lengths . 

INVE  "1 

FLOW 

TYPES 

I 

Numeric  code  for  Flow  key¬ 
word  =  16. 

FLPAR 

10,10 

CABLE 

R 

Flow  field  table  parameters 

FLOW  1 

FL'JI 

TYPES 

I 

Numeric  code  for  FLUI  key¬ 
word  =  17. 

FP 

3,NN3 

ACOM 

R 

Point  loads  for  load  sets. 

LOAD 

FRCFAC 

FRQDAT 

R 

Force  conversion  factor  for 
ship  motion  file. 

EXTE  1 

FRCVY 

CONTRL 

R 

Inventory  conversion  factor 
for  weights  and  strengths. 

INVE 

FREE 

TYPES 

I 

Numeric  code  for  FREE  key¬ 
word  =  30. 

FREQ 

TYPES 

I 

Numeric  code  for  FREQ  key¬ 
word  =  6. 

FSFRC 

3 

SHIPS 

K 

Ship  load  table  conversion 
factor  for  current  force. 

SHIP 

FSFRW 

3 

SHIPS 

R 

Ship  load  table  conversion 
for  wind  force. 

SHIP 

FSLEN 

3 

SHIPS 

R 

Ship  load  table  conversion 
factor  for  length. 

SHIP 

FSOL 

TYPES 

I 

Numeric  code  for  FSOL  key- 

i 

word  =  47. 
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VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYVORD 

FSVEL 

3 

SHIPS 

R 

Ship  load  table  conversion 
for  velocity. 

SHIP 

FTF 

3 

TIMED 

R 

Time  fimction^values  for 
load  sets  in  C. 

FTI 

15 

TIMED 

R 

Time  function  values  for 
dynamic  moving  boudary  com¬ 
ponents  in  C. 

FVISC 

5 

CABLE 

R 

Fluid  table  kinematic  vis¬ 
cosities. 

FLU  I 

FI 

NN3 

ACOM 

R 

Nodal  point  loads  frgi® 
current  and  wind  in  C. 

F2 

NN3 

ACOM 

R 

Total  nodal  point  loads  at 
the  incremental  reference 
state.  Also  used  for  itera¬ 
tive  displacement  increment 
in  NRAPIT  subroutine. 

GAIR 

SHIPS 

R 

Specific  weight  of  air. 

GAMNF.W 

TIMED 

R 

Gamma  integration  parameter. 

TIME 

GHED 

FRQDAT 

R 

Global  wave  heading  in 
degrees . 

HEAD 

GK 

NF3,IBEND 

ACOM 

R 

Global  incremental  stiff¬ 
ness  matrix  stored  in  banded 
row  form. 

GKS 

3,3 

SHIPS 

R 

Element  stiffness  matrix  for 
strum  string. 

GM 

NN3 

ACOM 

R 

Diagonal  mass  matrix  includ¬ 
ing  residual  added  mass. 

GMA 

NE 

ACOM 

R 

Residual  added  mass  for  each 
element. 

GKU 

6,6 

FRQDAT 

R 

Ship's  mass  matrix  from  ship 
motion  file. 

GRAV 

CONTRL 

R 

Acceleration  due  to  gravity. 

PR06 

Gl 

CONTRL 

R 

Kinematic  viscosity  of  first 
fluid,  water. 
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VARIABU: 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

G2 

CONTRL 

R 

Specific  weight  of  first 
fluid,  water. 

G3 

10 

CABLE 

R 

Weight  per  unit  length  of 
line  material. 

HATE 

HAWS I Z 

23 

COMPNT 

R 

Inventory  hawser  sizes. 

HAVTAB 

2,23,4 

COMPNT 

R 

Inventory  hawser  property 
table . 

HDGl 

30 

FRQDAT 

R 

Wave  headings  for  ship 
motion  file. 

HEAD 

3 

SHIPS 

R 

Ship's  global  heading. 

HED 

20 

HEDCHR 

H 

Page  heading  title. 

HEDEND 

SHIPS 

R 

Surface  load  total  heading 
change . 

STEP 

HEDliiX 

SHIPS 

R 

Upper  bound  on  heading 
increment. 

HEDINC 

SHIPS 

R 

Surface  load  heading 
increment. 

STEP 

HEDNOW 

SHIPS 

R 

Surface  load  heading  value 
in  ^C. 

HIRSDL 

DSPCON 

R 

VRR  control  parameter;  keeps 
tract  of  largest  residual 
found . 

HLDEAC 

6 

COMPNT 

R 

Inventory  anchor  holding 
factor  table. 

HTYPE 

3,4 

CWIPNT 

R 

Inventory  hawser  types. 

lAABU 

10 

IBUOYS 

I 

Additional  array  for  future 
expansion. 

lAACA 

10 

I CABLE 

I 

Additional  array  for  future 
expansion. 

lAACL 

10 

ICNTRL 

I 

Additional  array  for  future 
expansion. 
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VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

lAADS 

10 

lOSPCN 

I 

Additional  array  for  future 
expansion. 

lAAPO 

10 

I PAYOT 

I 

Additional  array  for  future 
expansion. 

lAASHP 

10 

ISHIPS 

I 

Additional  array  for  future 
expansion. 

lAASTM 

10 

I STRUM 

I 

Additional  array  for  future 
expansion. 

lAATIM 

10 

ITIMED 

I 

Additional  array  for  future 
expansion. 

lALTR 

ITIMED 

I 

Alternating  acceleration 
increment  flag  for  DYN  SAO. 

IB 

5 

ITIMED 

I 

List  of  dynamic  moving 
boundary  nodes. 

I RKND 

ICNTRL 

I 

Half  bandwidth. 

IBKNDl 

ICNTRL 

I 

Half  bandwidth  with  incre¬ 
mental  moving  boundary. 

lBEiro2 

ICNTRL 

I 

Half  bandwidth  without 
incremental  moving  boundary. 

IBU 

ICNTRL 

I 

Debug  output  flag 

OUTP 

TBGF 

FRQUAT , 

1 

Debug  otilput  flag  for 

FRQ  S.AO. 

RESU 

i.'S 

50 

I  BUOYS 

1 

Time  variation  codes  for 
surface  buoys. 

SBUO 

13'l! 

50 

I BUOYS 

I 

Nodes  where  bodies  are 
located . 

BLOC 

ii;>iciiF 

FRQDAT 

1 

Component  check  file  flag. 

FSOL 

J  CONMS 

FRQDAT 

I 

Mass  matrix  format  code. 

FSOL 

ic:k 


ISHIPS  I 


Flag  to  signal  if  table 
lookup  is  used  to  get  hull 
resistance  factor. 


VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

IDAHPR 

ICABLE 

1 

Material  damping  flag; 

0  means  no  damping. 

lOAY 

HEDCHR 

H 

Character  string  for  date. 

IDIR 

ICNTRL 

1 

Coordinate  number  for  verti¬ 
cal  direction. 

IDIS 

30,3 

IDSPCN 

I 

Static  displaceotent  variation 
codes . 

IDRB 

50 

I BUOYS 

I 

Drag  coefficient  codes  for 
bodies . 

BODY 

IDRG 

10 

ICABLE 

I 

Drag  coefficient  codes  for 
line  materials. 

MATE 

IDRITR 

FRQDAT 

I 

Drift  force  iteration 
flag. 

FSOL 

IDWN 

IDSPCN 

I 

VRS/VRR  control  parameters; 
signals  successive  reduc¬ 
tions  of  velocity  norm. 

IFCNT 

ICNTRL 

I 

Record  counter  for  frequency 
domain  CHEK  data. 

IFILE 

4 

SIZE 

I 

Restart  file  record  counters. 

IFLAGS 

60 

IFLAG 

H 

Codes  for  each  of  the  key¬ 
words,  4  characters  each. 

IFLCOD 

10 

ICABLE 

I 

Flow  field  codes. 

FLOW 

IFREQ 

I SHIPS 

I 

Signals  the  frequency  domain 
solution  to  the  stiffness 
matrix  routines. 

IFRQUP 

FRQDAT 

I 

Iteration  option  for  drift 
force  updates. 

FSOL 

IFXFL 

ICNTRL 

I 

Limit  condition  flag  used 
to  signal  constraint  code 
changes . 

ICK 

ICNTRL 

I 

Solution  format  code  for 
stiffness  matrix. 
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VARIABLE 

NAME 

ARRAY 

sm: 

COMMON 
NAME  . 

TYPE 

DESCRimON 

INPUT 

KEYWORD 

I  KEEP 

TYPES 

I 

Numeric  code  for  KF,KP  key¬ 
word  =  33. 

IKNSTN 

ICNTRL 

I 

Limit  check  over-ride  flag. 

ILF 

i 

ITIMED 

I 

Load  set  variation  codes. 

LVAR 

IMPBOD 

I BUOYS 

I 

Body  mwnher  for  body  impact. 

IMPA 

IMPNOD 

I BUOYS 

I 

Node  number  for  body  impact. 

IMPA 

IMTMF 

15 

ITIMED 

I 

Moving  boundary  variation 
codes . 

MOVE 

IMX 

ITIMED 

I 

Number  degree  of  freedom 
with  largest  acceleration 
increment  for  DYN  SAO. 

INC 

ICNTRL 

1 

Solution  increaient  counter. 

INDRAG 

ICNTRL 

I 

Drag  model  override  flag. 

PROB 

INVFLG 

CHKDAT 

I 

Inventory  useage  flag. 

INVP 

IDSPCN 

I 

VRS/VRR  control  parameter. 

TNW 

ICNTRL 

I 

Inventory  useage  flag. 

lOPT 

I BUOYS 

1 

Impacting  body  weight  option 
code. 

IMPA 

i  r’AGE 

IIEDDAT 

I 

Page  counter. 

Il'R 

ICNTRL 

1 

Number  of  steps  between 
printouts . 

OUTP 

j  KOI  IT 

FRQDAT 

I 

Roll  damping  iteration  flag. 

FSOL 

IRK 

I 

Error  return  flag;  means 

error. 

IRST 

RETAPE 

I 

Restart  file  flag. 

SAVE 

or 

REST 

ISHIP 

:i 

I SHI PS 

I 

Nodes  where  ships  are 

SHIP 

located . 
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VAKIABLK 

NAME 

ARRAY 

SIZE 

COMMON 

namf: 

TYPE 

DESCRIPTION 

iNPirr 

KEYWORD 

ISHPFL 

FRQDAT 

1 

Ship  flag  for  FREQ  SAO; 

0  means  no  ship. 

ISHTAP 

[SHIPS 

1 

Flag  to  signal  the  use  of 
ship  load  file. 

ISIGNI 

ITIMED 

1 

Alternating  response  esti¬ 
mates  flag  for  DYN  SAO. 

ISIGN2 

[TIMED 

I 

Alternating  response  esti¬ 
mates  flag  for  DYN  SAO. 

1ST 

CDCAI, 

I 

Strum  string  niunbcr. 

I START 

[CNTRL 

I 

Solution  option  code;  see 
NSLOP. 

ISTRNG 

30 

1 STRUM 

I 

Number  of  elements  in  each 
strum  string. 

STRU 

ISTRUP 

I STRUM 

I 

St  rum  update  flag. 

ISURLD 

ISHIPS 

I 

Wind  and  surface  current 
load  flag. 

IT 

2,NK 

ACOM 

1 

Element  connectivity  list 

ELEM 

ITCONF 

FRQDAT 

I 

Counter  for  drift  force 
iterations , 

ITFCOI) 

20 

ITIMED 

I 

Time  function  type  code. 

TFIIN 

rriME 

HEDCHR 

H 

Character  string  for  time. 

ITOP 

ITIMED 

1 

Upper  bound  flag  on  time 
step  in  DYN  SAO. 

ITIIF 

inSPCN 

I 

VRS/VRR  solution  parameter. 

IDNKES 

FRQDAT 

I 

Unrestrained  motion  ou..put 
flag. 

KESU 

11)1‘ 

ISHIPS 

I 

Code  far  up  direction. 

1 UPOT 

1 CNTRL 

I 

Counter  for  configuration 
updates. 

JANCB 

*>0 

I BUOYS 

1 

Limit  set  fixity  codes. 

L.TMI 
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VARIABl.K 

NAMK 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

JB 

15 

1 TIMED 

1 

Moved  component  variation 
codes . 

MOVE 

JDt.D 

ICNTRL 

1 

Dead  load  option  flag 

LVAK 

JDYN 

ICNTKL 

I 

Dynamic  solution  option 
flag. 

PROB 

JMPDT 

ICNTRL 

I 

Step  size  control  niunbcr. 

SOLU 

JOP 

5 

1  PAYOT 

I 

Node  numbers  for  payout 
points. 

PAYO 

JOVR 

ICNTRL 

1 

Limit  set  overshoot  node. 

JPEl.T 

5 

U’AYOT 

1 

Element  numbers  for  payout  . 

PAYO 

JSLP 

2,50 

1  BUOYS 

I 

Lists  connecting  elements 
for  cylindrical  buoys  or 
defining  nodes  for  mooring 
buoys . 

JSTEPR 

IDSPCN 

1 

Niunber  of  steps. 

STEP 

KNSTHN 

IBUOYS 

I 

Flag  to  signal  if  limit 
locations  arc  defined. 

KN'^P'^T 

ITTMED 

I 

Counter  for  restart  save 
i nterval s . 

KODtC 

•) 

1  SHI  PS 

1 

Variation  codes  for  surlace 
current . 

fiURF 

KODEv, 

7 

1 SHIPS 

I 

Variation  codes  for  wind 
loads . 

WIND 

NE 

ACOM 

I 

Element  type  code. 

ELEM 

KONKi T 

10,50 

IHUOYS 

I 

Lists  elements  connecting 
to  each  limit  location. 

KiAVX'l' 

TCNTRI, 

1 

Number  of  iteration  trials. 

SOLU 

KOUNT 

ICNTRL 

I 

Solution  iteration  counter. 

KSTRNG 

20,30 

I  STRUM 

1 

Lists  strum  string  elements 
up  to  20  per  string. 

STRU 
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VARIABLE 

NAHF. 

ARRAY 

COMMON 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

KLIP 

ICNTRL 

1 

Time  step  increase  control 
parameter  in  DYN  SAG. 

LBODN 

50 

1 BUOYS 

I 

List  of  body  numbers  of 
body  locations. 

BLOC 

LIMNOl) 

50 

1 BUOYS 

I 

List  of  limit  location 
nodes . 

U.0C 

LIMSET 

50 

I BUOYS 

1 

List  of  limit  set  numbers 
for  limit  locations. 

LI.0C 

LLBI. 

SHPLBL 

H 

Length  label  for  ship  load 
file. 

RIGID 

LiiiTKR 

ICNTRL 

I 

Maximiun  number  of  itera¬ 
tions  per  step. 

SOLU 

I^KEEP 

IDSFCN 

I 

Temporary  storage  of 

LMITER  for  VRS/VRR 
solutious . 

LSHl> 

3 

1 SHI PS 

1 

Ship  load  function  codes. 

SHIP 

MAT 

NE 

ACOM 

1 

List  of  elcmeut  inalerial 
number.s. 

ELEM 

MATDMl' 

10 

1 CABLE 

I 

Li.st  of  material  damping 
flags;  0  means  this  material 
has  no  damping. 

MATT 

1 CABLE 

1 

Number  of  materials  defined. 

MBLN 

50 

T BUOYS 

I 

Counters  for  number  of 
number  slave  nodes  on  aioor- 
ing  buoys. 

MK 

10 

I CABLE 

1 

Number  of  points  in  each 
material  table. 

MATE 

MED 

10 

1 CABLE 

1 

List  of  medium  codes  for 
material  tables. 

MATE 

MED  HIM 

NF. 

ACOM 

I 

Fluid  mediiun  codes  for 
elements . 

ELKM 

MEOMB 

50 

I BUOYS 

I 

List  vif  medium  <'<>d«'s  for 
body  table. 

BODY 
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VARIABLE  ARRAY 
NAME  SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYVORO 

MITNOT 

I PAYOT 

I 

Mitosis  preventioa  flag. 

MNOB 

FRQDAT 

I 

Maxiffluoi  niunber  of  wave 
speeds  on  ship  notion  file. 

MNOH 

FRQDAT 

I 

Maximum  number  of  wave  head* 
ings  on  ship  motion  file. 

MNOK 

FRQDAT 

I 

Maximum  number  of  wave 
lengths  on  ship  motion  file. 

MNRV 

FRQDAT 

1 

Maximum  number  of  roll 
angles  on  ship  motion  file. 

MODE 11 

ICNTRL 

1 

Mode  shape  order  flag. 

MSOL 

MODE I 2 

ICNTRL 

I 

Mode  shape  output  flag. 

MSOL 

MORBUO 

I BUOYS 

I 

Number  of  mooring  buoys. 

MOSTAT 

ICNTRL 

I 

TSSS  flag;  1  -  yes,  o  -  no. 

MULMAT  5 

i PAYOT 

I 

Multimaterial  flag  for 
payout  element. 

MULTIM 

I PAYOT 

I 

Multimaterial  payout  flag. 

MVB 

ICNTRL 

I 

Number  of  nodes  with  defined 
motion  in  DYN  SAG. 

MV B INC 

ICNTRL 

I 

Number  of  incremental  moving 
boundary  nodes. 

KXBLOC 

I BUOYS 

I 

Maximum  number  of  body  loca¬ 
tions. 

MXBODY 

I BUOYS 

1 

Maximum  number  of  bodies 
in  table. 

MXFLOW 

1 CABLE 

I 

Maximum  number  of  flow 
tables. 

MXLIMS 

IBUOYS 

I 

Maximum  number  of  limit 
sets. 

MXLINE 

HEDDAT 

I 

Number  of  lines  per  page. 

MXLLOC 

IBUOYS 

I 

Maximum  number  of  limit 
locations. 
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VARIABLE 

ARRAY 

COMMON 

NAME 

SIZE 

NAME 

TYPE 

MXSELT 

ISTRUM 

I 

HXSUIP 

ISHIPS 

I 

MXSTRG 

ISTRUM 

I 

MXTFUN 

ITIMED 

I 

NANCR 

6 

COMPNT 

I 

NBASE 

ICNTRL 

I 

NBLOC 

IBUOYS 

I 

NBUOY 

2 

COMPNT 

I 

NCHN 

COMPNT 

I 

NCOM 

SIZE 

I 

NCONC 

ICNTRL 

I 

NCRNT 

ISHIPS 

I 

NCYLB 

IBUOYS 

I 

NDISP 

IDSPCN 

I 

NDROP 

IDSPCN 

1 

N£ 

ICNTRL 

I 

NELCK 

CHKDAT 

1 

NELPOl 

5 

I PAYOT 

I 

INPUT 

DESCRIPTION  KEYWORD 

Maximum  number  of  element 
per  gtruB  string. 

Maximum  number  of  ship 
locations. 

Maximum  number  of  strum 
strings . 

Maximum  number  of  time 
functions. 

Inventory  anchor  count. 

Number  of  entries  in  the 
base  portion  of  ACOM. 

Number  of  body  locations. 

Inventory  buoy  count. 

Inventory  chain  count. 

Size  of  ACOM. 

Number  of  concentrated  loads 
(point  loads). 

Number  of  current  tables  in 
ship  load  file. 

Number  of  cylindrical  buoys. 

Niusber  of  static  imposed 
displacement. 

VRS/VRR  control  parameter; 
counts  down  for  3  velocity 
convergence  trials  with 
damping  reduced. 

Number  of  cable/line  ele¬ 
ments  . 

Element  number  being 
checked . 

Element  number  increment  PAYO 

for  payout. 
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VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

NEWCON 

FRQDAT 

I 

New  configuration  flag 
for  drift  force  updates. 

NEWRED 

FRQDAT 

I 

Global  heading  input  flag 
for  FREQ  SAO. 

NFILE 

RETAPE 

I 

Restart  file  code 

REST 

NFILEF 

FRQDAT 

I 

Ship  motion  file  format 
flag. 

EXTE 

NFIX 

ICNTRL 

1 

Number  of  globally  fixed 
nodes  (always  0) . 

NFKEEP 

30 

IDSPCN 

I 

Previous  nodes  component 
fixity  codes  for  static 
imposed  displacement  and 
moving  boundary. 

NFLUl 

I CABLE 

I 

Number  of  fluids  defined. 

NFI.VRY 

I CABLE 

I 

Flow  field  variation  code 

CURR 

KFLUID 

ICNTRL 

I 

Flow  field  number 

CURR 

NFN 

ICNTRL 

I 

Number  of  active  nodes. 

ICNTRL 

I 

Number  of  active  degrees 
of  freedom. 

is'FKEQS 

FRQDAT 

1 

Counter  for  the  number  of 
frequencies  in  FREQ  SAO. 

N(}ROW 

!) 

I PAYOT 

1 

Number  of  elements  available 
for  payout . 

PAYO 

N'HAWS 

COMPNT 

I 

Inventory  hawser  count . 

NliED 

liEDDAT 

I 

Number  of  words  in  page 
heading;  defined  in  SEADYN. 

NIN 

TAPES 

I 

Input  file  code. 

NINA 

SIZE 

1 

Number  of  words  in  /ACOM/. 

N  INB 

SIZE 

I 

Number  of  words  in  /BUOYS/. 

NINRI 

SIZE 

I 

Number  of  words  in  /IBUOYS/. 
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VAKIAHLK 

ARRAY 

.yZE 

COMMON 

JIAfl|K_ 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

NINC 

SIZE 

1 

Number  of  words  in  /CABLE/. 

NINCI 

SIZE 

I 

Number  of  words  in  /ICABLE/. 

NlNCLl 

SIZE 

1 

Number  of  words  in  /ICNTRL/. 

NlNBSl 

SIZE 

I 

Niunber  of  words  in  /IDSPCN/. 

NINDSP 

SIZE 

1 

Number  of  words  in  /DSPCO.V/. 

Nl.NPO 

SIZE 

I 

Number  of  word.s  in  /PAYOUT/. 

NINPOI 

SIZE 

I 

Number  of  words  in  /IPAYOT/. 

NINP.L 

SIZE 

I 

Number  of  words  in  /CONTRL/. 

NlNSHl 

SIZE 

I 

Number  of  words  in  /ISHIPS/. 

NINSHP 

SIZE 

I 

Number  of  words  in  /SHIPS/. 

N INST  I 

SIZE 

T 

Number  of  words  in  /I STRUM/. 

NTNSTM 

SIZE 

1 

Number  of  words  in  /STRUM/. 

NINT 

SIZE 

1 

Number  of  words  in  /TIMED/. 

NlNT.Ml 

SIZE 

I 

Number  of  words  in  /ITIMED/. 

N 1 XPPN 

/TIMED 

1 

Print  override  flaji;  for 
restart . 

SAVE 

Ni.;;is 

I  BUOYS 

I 

Number  of  limit  sets. 

NLiNES 

llEDDAT 

I 

Output  lines  counter. 

Nl-LOC 

1  BUOYS 

I 

Number  of  limit  locations. 

NMOTN 

.SO 

I  BUOYS 

1 

Buoy  on  surface  motion  code. 

SBUO 

MN 

ICNTRL 

1 

Total  number  of  node.s. 

NNPDl 

S 

I  PAYOT 

I 

Not  used  at  present 

NN;) 

ICNTRL 

I 

Three  times  NN. 

NOB 

FRQDAT 

1 

Number  of  wave  speeds  on 

ship  motion  file. 
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VARIABLE 

NAME 

NODFIX 

NODWN 

NOFLUD 

NOH 

NOITER 

NOK 

NOLINE 

NOLOAD 

NOP 

NOUT 

NOVEL 

NPOVRY 

NPRECZ 

NPKST 

NRUP 

NRV 

NSFILE 

NShIPS 

NSHRNK 


ARRAY 

SIZE 

NN3 


5 


5 


COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

ACOM 

I 

Node  component  fixity  flags. 

IDSPCN 

I 

VRS/VRR  control  parameter. 

LOGIC 

L 

Fluid  medium  flag;  True 
means  no  fluids . 

FRQDAT 

I 

Number  of  wave  headings  on 
ship  motion  file. 

LOGIC 

L 

Iteration  flag;  T  means  not 
iterative  solution. 

FRQDAT 

I 

Number  of  wavelengths  on 
ship  motion  file. 

HEDDAT 

I 

Number  of  lines  is  standard 
output  record. 

LOGIC 

L 

Point  load  flag;  True  means 
no  point  loads. 

IPAYOT 

I 

Number  of  payout  ends. 

TAPES 

I 

Output  file  number. 

LOGIC 

L 

Flow  field  flag;  True  means 
no  fluid  velocities. 

IPAYOT 

I 

Payout  time  variation  codes. 

PAYO 

SIZE 

1 

Precision  of  floating  point 
words;  1  -  single,  2  -  double. 

ICNTRL 

I 

Preload  flag. 

ICNTRL 

I 

Newton-Raphson  update  code. 

SOLU 

FRQDAT 

I 

Number  of  roll  angles  on 
ship  motion  file. 

ISHIPS 

I 

Ship  load  file  flag. 

PROS 

ISHIPS 

I 

Number  of  ships. 

IPAYOT 

I 

Number  of  elements  avail¬ 
able  for  reel-in 

PAYO 
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VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT  13 

KEYWORD 

NSLAVE 

ICNTRL 

I 

Nuttber  of  slave  nodes. 

1:‘ 

Vi 

1( 

NSLOP 

IDSPCN 

I 

Solution  option  number  code; 

0  -  MNR,  1  -  SLI,  -  RFB, 

3  -  VRS/VRR,  5  -  DIM. 

-  -I.i 

NSOLN 

FRQDAT 

I 

Counter  for  FREQ  wave  head¬ 
ing  solutions. 

J 

NSRDC 

FRQDAT 

I 

File  code  for  ship  motion 
file  =  8. 

NSTRNG 

I STRUM 

I 

Number  of  strum  strings. 

NSTRUP 

I STRUM 

I 

Strum  update  flag. 

NSTUP 

ICNTRL 

I 

First  step  subdivision 
startup  parameter. 

STEP 

NS3 

FRQDAT 

I 

Three  times  (NN-NSLAVF.S) . 

\ 

; 

NTAPE 

RETAPE 

I 

File  code  for  restart  input 
tape . 

REST 

NTAPEl 

TAPES 

I 

File  code  for  binary  trans¬ 
lation  of  free-field  input 
-  15. 

NTAPE2 

TAPES 

1 

File  code  for  rigid  format 
input  =  16. 

1 

NTAPE9 

TAPES 

I 

File  code  for  file  code  09, 
scratch  file. 

:■  ^ 

NTHETC 

ISHIPS 

I 

Number  of  current  headings 
in  ship  load  tables. 

RIGID  1 

NTllETW 

I  SHIPS 

I 

Number  of  wind  headings 
in  ship  load  tables. 

RIGID 

\ 

NTPCHK 

TAPES 

I 

File  code  for  wave  heading 
file  for  CHEK  -  13. 

1 

NTYPE 

ICNTRL 

I 

Numeric  code  for  current 

SAO  keyword. 

NUmiED 

ISHIPS 

I 

Counter  for  the  number  of 

headin^K  evaluated  in  static 
heading  excursion. 
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VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

NUMSET 

ICNTRL 

I 

Nunber  of  load  sets  defined. 

NUP 

ICNTRL 

I 

Configuration  update  flag 

SOLU 

NWINO 

SHIPS 

I 

Number  of  wind  velocity 
tables  on  ship  load  file. 

RIGID 

NWORDS 

HEDDAT 

I 

Number  of  words  in  working 
area  for  input. 

NXTYPE 

ICNTRL 

1 

Numberic  code  for  the  next 

SAO  keyword. 

OMG 

CDCAL 

R 

Strum  string  natural 
frequency. 

OMG 

FRQDAT 

R 

Spectrum  frequency  for 

FREQ  SAO. 

OMGHN 

FRQDAT 

R 

Lower  bound  of  spectrum 
frequency  .scan 

SPEC 

OMGMX 

FRQDAT 

R 

Upper  bound  of  spectrum 
frequency  scan. 

SPEC 

CMS 

CDCAL 

R 

Strum  string  Strouhal 
frequency. 

0M2 

FRQDAT 

K 

Square  of  OMG  for  FRQ  SAO. 

OIJTP 

TYPES 

I 

Numeric  code  for  OUTP  key¬ 
word  =  38. 

CVKSHT 

CONTRL 

R 

Limit  condition  overshoot. 

PARMT 

CONTRL 

R 

MNR  extrapolation  parameter. 

SOU) 

PAVO 

TYPES 

I 

Numeric  code  for  PAYO  key¬ 
word  =  39. 

PAYV 

5 

PAYOUT 

R 

Payout  velocity 

PAYO 

PI 

CONTRL 

R 

PI  =  3.1415 . 

PINC 

CONTRL 

R 

Output  interval  increment. 

PLUT 

TYPES 

1 

Numeric  code  for  PLOT  key¬ 
word  =  8,  inactive. 
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VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

POUT 

LOGIC 

L 

Payout  flag;  True  neaos 
payout  active. 

PROB 

TYPES 

I 

Numeric  code  for  PROB  Key** 
word  =11. 

PROPF 

3 

SHIPS 

R 

Propeller  force  coefficients. 

PSTEP 

CONTRL 

R 

Output  load  level. 

PTMF 

5 

PAYOUT 

R 

Values  of  payout  time 
variation  functions . 

QST 

20 

STRUM 

R 

Strum  string  mode  shape. 

RAND 

TYPES 

I 

Numeric  code  for  RAND  key¬ 
word  =  52. 

RANG 

8 

FRQDAT 

Roll  angles  for  ship  motion 
file. 

RAID 

CONTRL 

R 

Initial  size  of  static  load 
increment. 

RAIL 

3 

SHIPS 

R 

Length  ratios  for  ship's 
current  load  scaling. 

RRFUP 

LOGIC 

L 

Geometric  reference  update 
flag. 

REGU 

TYPES 

I 

Numeric  code  for  REGU  key¬ 
word  =  53. 

RELFAC 

50 

BUOYS 

R 

Release  factor  for  limit 
sets. 

RKRR 

CONTRL 

R 

Residual  error  tolerance. 

REST 

TYPES 

I 

Numeric  code  for  REST  key¬ 
word  =  12. 

RESU 

TYPES 

1 

Numeric  code  for  RESU  key¬ 
word  =  50. 

RNORM 

DSPCON 

R 

Residual  norm. 

INPUT 

KEYVORD 


PAYO 


UMI 

SOLU 


RNORMP 


OSPCON  R 


Preceding  residual  nom. 


VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

RNRMPP 

DSPCON 

R 

Second  preceding  residual 
norm. 

RTIN 

CDCAL 

R 

Strum  string  modal  scaling 
factor. 

RVELN 

STRUM 

R 

Relative  velocity  magnitude 
for  strum  computation. 

RVELNP 

STRUM 

R 

Previous  relative  velocity 
magnitude  for  strum  compu¬ 
tation. 

SAE 

3 

SHIPS 

R 

Ship  end  projected  areas. 

SHIP 

SAS 

3 

SHIPS 

R 

Ship  side  project  areas. 

SHIP 

SAVE 

TYPES 

I 

Numeric  code  for  SAVE  key¬ 
word  =  40. 

SBAMP 

50 

BUOYS 

R 

Surface  buoy  motion  ampli¬ 
tudes. 

SBUO 

SBEAM 

3 

SHIPS 

R 

Ship  beams  amidships. 

SHIP 

SBUO 

TYPES 

1 

Numeric  code  tor  SBUO  key¬ 
word  =  4l. 

SCALE 

SHPTAP 

R 

Scale  factor  for  ship  load 
file. 

RIGID 

SDRFT 

3 

SHIPS 

R 

Ship  drafts  amidships. 

SHIP 

SDSPV 

> 

SHIPS 

R 

Ship  volume  displacements. 

SHIP 

SfACW 

3,3 

SHIPS 

R 

Wind  load  scale  factors 
for  ships . 

Si.NDLD 

BUOYS 

R 

Gravity  load  sign;  +1.0  for 
gravity  in  +  coordinate  di¬ 
rection,  -1.0  for  gravity 
in  -  coordinate  direction. 

vSHPCAP 

12 

SHPLBL 

H 

Title  for  ship  load  file. 

RIGID 

Ship's  hydrostatic  restor¬ 
ing  coefficients. 


SHTPK 


4,3 


SHIPS 


K 


SHIP 


VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPJ 

DESCRIPTION 

INPUT 

KEYWORD 

SHPKP 

3,3 

SHIPS 

R 

Ship  response  parameters 
for  checking  static  response 
stability. 

SHTRN 

3,3,3 

SHIPS 

R 

Local  to  global  transfocna- 
tion  matrices  for  each  ship. 

SIG 

N£ 

ACOM 

R 

Element  tensions  in  ^C. 

ELEM  or 
TENS 

SIGR 

NE 

ACOM 

R 

D 

Element  tensions  in  C. 

SLT 

3 

SHIPS 

R 

Ship  lengths. 

SHIP 

SLWL 

3 

SHIPS 

R 

Ship  water  line  lengths. 

SHIP 

SOI.U 

TYPES 

I 

Numeric  code  for  SOLU  key¬ 
word  ®  A2 . 

SPEC 

TYPES 

I 

Numeric  code  for  SPEC  key¬ 
word  =  48. 

SPECA 

FRQDAT 

R 

Wave  spectrum  coefficient,  A. 

SPEC 

SPECb 

FRQDAT 

R 

Wave  spectrum  coefficient,  B. 

SPEC 

SRCHKC 

CONTRL 

R 

Search  factor  for  II)  search 
accelerator  on  MNR  solution. 

SOLU 

SSTART 

CONTRL 

R 

First  step  subdivision 
parameter. 

STEPUP 

LOGIC 

L 

Increment  reference  update 
flag. 

STLEN 

20 

STRUM 

R 

Strum  string  lengths. 

STN 

N£ 

ACOM 

R 

Element  strains  in  ^C. 

st:: 

20,10 

CABIJ; 

R 

Strain  values  for  material 
tables . 

MATE 

SURF 

TYPES 

I 

Numeric  code  for  SURF  key¬ 
word  =  44. 

SURFCE 

SHIPS 

R 

Vertical  coordinate  at  the 

water  surface,  detertniued 
by  the  initial  vertical  posi¬ 
tion  of  the  last  ship  input. 
Used  only  when  ship  restor¬ 
ing  coefficients  are  input. 
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VARIABLE 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

T 

TIMED 

R 

Time. 

TBLOCK 

SHPTAP 

R 

Test  Ship  block  coefficient 
for  ship  load  file. 

TDEPTH 

SHPTAP 

R 

Test  depth  for  ship  load 
file. 

TEMP 

6,6 

F8QDAT 

R 

Temporary  storage  array  for 
FREQ  SAO. 

TENS 

TYPES 

1 

Numeric  code  for  TENS  key¬ 
word  “  26. 

TENULT 

SO 

CALBE 

R 

Line  material  ultimate 
tensions. 

MATE 

TFSAV 

5 

PAYOUT 

R 

Temporary  storage  of  payout 
time  functions. 

TFUN 

TYPES 

* 

Numeric  code  for  TFUS  key¬ 
word  =  27 . 

TN 

3,NE 

ACOH 

R 

Direction  cosines  for  each 
element  in  C. 

THR 

3,NE 

ACOM 

R 

Direction  cosines  for  each 
t  .eat  in  C. 

TIMFAC 

FRQDAT 

R 

Time  conversion  factor  for 
ship  motion  file. 

EXTE 

TMAS 

FRQDAT 

R 

Ship's  mass  for  ship  motion 
file. 

TMAX 

TIMED 

R 

Maximum  time. 

TMFRF 

3 

TIMED 

R 

Point  load  tige  variation 
functions  in  l. 

TMFRM 

15 

TIMED 

R 

Moving  boundary  time  » 
variation  functions  in  c. 

miM 

SO 

BUOYS 

R 

Limit  set  tolerances. 

LIMI 

TPARM 

20,20 

TIMED 

R 

Tin;e  function  tables. 

TFUN 

TRNSTR 

3,3,NE 

ACOH 

R 

l.ocal-to-global  transforma¬ 
tion  matrices  for  each 
element. 
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VAKIABU; 

NAME 

ARRAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INHJT 

KEYWORD 

TKSRT 

TIMED 

R 

Restart  save  time. 

TSAK 

SHPTAP 

R 

Test  ship  end  projected 
area  for  ship  load  file. 

RIGID 

TSAP 

SHPTAP 

R 

Test  ship  propeller  pro¬ 
jected  area  for  ship  load 
file. 

RIGID 

TSAS 

SHPTAP 

K 

Test  ship  size  projected 
area  for  ship  load  file. 

RIGID 

TSAPL 

3 

SHIPS 

R 

Test  ship  propeller  pro¬ 
jected  areas;  used  for 
similarity  scaling. 

TSB 

SHPTAP 

K 

Test  ship  beam  for  ship 
load  file. 

RIGID 

TSl) 

SHPTAP 

R 

Test  ship  draft  for  ship 
load  file. 

RIGID 

TSDSP 

SHPTAP 

R 

Test  ship  displacement  for 
ship  load  file. 

RIGID 

TSLT 

SHPTAP 

R 

Test  ship  length  for  ship 
load  file. 

RIGID 

TSSS 

TYPES 

I 

Numeric  code  for  TSSS  key¬ 
word  =  A. 

TSVL 

SHPTAP 

R 

Test  ship  waterline  length 
for  ship  load  file. 

RIGID 

■IT 

7.0,10 

CABl.E 

R 

Teusion  values  for  material 
tables . 

MATE 

TTD 

10 

CABLE 

R 

Material  damping  parameter 

CA^. 

MATE 

TTK 

10 

CABLE 

R 

Material  stiffness  parameter 
EA,. 

MATE 

TIMED 

K 

Time  when  incremental  refer¬ 
ence  was  established. 

11 

NN3 

ACOM 

R 

Nodal  displacement  incre- 

uonts  irom  in 
reference  to 
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VAKIABLF. 

NAME 

AKKAY 

SIZE 

COMMON 

NAME 

TYPE 

DESCRIPTION 

INPUT 

KEYWORD 

l)B 

IS 

TIMED 

H 

Moving  boundary  motion 
amplitudes . 

MOVE 

UBS 

30 

BUOYS 

R 

Surface  buoy  motion 
amplitudes. 

SBUO 

l)D 

NN3 

ACOM 

K 

Nodal  point  velocities 
for  ’•C. 

HDD 

NN3 

ACOM 

R 

Nodal  point  acclerations 
for  T,. 

UDDP 

NN3 

ACOM 

R 

Noda|_g^iut  accelerations 
for  C. 

UDDS 

NNE 

ACOM 

K 

Nodal  point  accelerations 
lor  incremental  reference 
state . 

UUP 

NN3 

ACOM 

R 

point  velocities  for 

X/  * 

UMVR 

IS,  3 

PAYOUT 

R 

Moving  boundary  motion  data 
for  '^C. 

lIMVBP 

IS,  3 

PAYOUT 

K 

Movijpg^^oundary  motion  data 
for  . 

UP 

NN3 

ACOM 

R 

Noda^_j^pnt  displacements 
for  C. 

US 

NN3 

ACOM 

R 

Nodalj.point  displacements 

1 rora  "C  to  incremental 
reference  state. 

VAIR 

SHIPS 

R 

Kinematic  viscosity  of  air. 

VDUM 

3 

FKQUAT 

R 

Temporary  storage  vector 
for  FREQ  SAG. 

VF 

NN3 

ACOM 

R 

Nodal  point  compo^ients  of 
flow  velocity  in  C. 

VI  B 

3 

BUOYS 

R 

Components  of  impacting 
body  velocity. 

I  MPA 

VI, Bl, 

SHPLBL 

H 

Velocity  label  for  ship 
load  file. 

RIGID 
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VARIABLE 

NAME 

ARRAY 

COMHUN 

NAME 

TYPE 

DESCRIPTION 

INPUT 

VNORMP 

DSPCON 

R 

Preceding  velocity  norm 
tor  VRS/VRR. 

VNPP 

DSPCON 

R 

Second  preceding  velocity 
norm. 

VNPPP 

DSPCON 

K 

Third  preceding  velocity 
norm. 

V.i 

NN3 

ACOM 

R 

Nodal  point  components  of 
relative  velocity  in  ‘C. 

WAD 

SHIPS 

R 

Wind  heading  in  degrees 

WIND 

WAVEL 

30 

FRQDAT 

R 

Wavelengths  from  ship 
motion  file. 

WDEPTH 

3 

SHIPS 

K 

Water  depth  for  each  ship. 

SHIP 

WIND 

SHIPS 

R 

Wind  velocity 

WIND 

WLBL 

SHPLBL 

H 

Wind  label  for  ship  load 
file. 

RIGID 

Wl.D 

LOGIC 

1, 

Live  load  flag;  True  means 
LIVE  SAO. 

WINDCOE 

20,3,5 

SHIPS 

R 

Wind  load  coefficients  for 
ship  load  table. 

RIGID 

WNDHKD 

20 

SHIPS 

H 

Wind  headings  for  ship 
load  tables. 

RIGID 

wNDVEL 

5 

SHIPS 

R 

Wind  velocities  for  ship 
load  tables. 

RIGID 

WORDS 

100 

HEDDAT 

R 

Working  area  for  input. 

WVAMP 

mQDAT 

R 

Wave  amplitude  for  current 
spectrum  interval. 

WVLN 

FRQDAT 

K 

Wave  length  for  current 
spectrum  iviterval. 

WVSLP 

P'RQDAT 

R 

Wave  .slope  for  current 
spectrum  interval . 

XC 

NN3 

ACOM 

R 

Nodal  position  coordinates 
for  C. 
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VARIABLE 

ARRAY 

COMMON 

NAME 

SIZE 

NAME 

TYPE 

XO 

>'N3 

ACOM 

R 

YSAVE 

CDCAL 

R 

INPUT 

DESCRIPTION  KEYWORD 

Reference  nodal  coordinates. 

Strum  string  response,  Y/D. 


9.0  MACRO-FLOW  CHARTS  OF  THE  SUBANALYSES  AND  SOLUTION  OPTIONS 


The  following  charts  outline  the  o'^erall  logic  of  the  major  solution 
options  in  SEADYN.  The  relationship  of  the  subanalyses  is  shown  in 
figure  2.1.  The  DEAD.  LIVE,  and  TSSS  subanalyses  are  controlled  from 
the  STATIC  subroutine.  Figure  9.1  outlines  that  logic.  The  transient 
dynamics  (DYN)  is  controlled  through  the  TRANS  subroutine.  That  logic 
is  pr;>sented  in  Figure  9.2.  Figure  9.3  describes  the  MODE  subanalysis. 
The  FREQ  and  CHEK  subanalyses  are  presented  in  Figures  9. A  and  9.5, 
respectively.  The  two  static  solution  routines  (NRAPIT  and  STEP)  are 
outlined  in  Figures  9.6  and  9.7. 
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Figure  9.3  Mode  subanalysis. 
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Figure  9-4.  Frequency  subanalytis. 
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Figure  9^.  Continued. 
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F|g*ic9.7  Stcplogie. 
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